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SYNOPSIS 

Drop size distributions in liquid-liquid dispersions 
evolve due to the dynamics of breakage and coalescence. 
Evolution of drop size distribution is important for 
prediction of, interfacial area, heat and/or mass transfer 
in liquid-liquid dispersions. In lean liquid -liquid 
dispersions, coalescence is negligible due to negligible 
interaction between droplets. The present investigation 
attempts to study the breakage of droplets in agitated lean 
liquid-liquid dispersions. 

Drop breakage in agitated lea,n liquid -liquid dispersions 
is random and hence is described by related probabilities, 
namely, the transitional breakage probability and the daughter 
droplet size distribution. The present investigation consists 
of three parts: 

(i) A model for transitional breakage probability 
(ii) A model for oquilibriim drop size distribution 
from a statistical thermodynamic viewpoint, 
and (iii) Experimental measurement of transient as well as 
equilibrium drop volume distributions in a batch 
agitated lean liquid-liquid dispersion in order 
to verify the two models and also to obtain the 
daughter droplet size distribution. 

The model for transitional breakage probability is 
based on the mechanism of breakage of droplets due to their 



oscillations resulting from relative velocity fluctuations, 

A universal transitional breakage probability in terms of 
non-dimensionalized drop volume is derived for all dispersed 
plaases whose density and viscosity are almost same as that 
of the continuous phase. The maximum stable drop diameter 
derived from this model shows a dependence of I'fc is 

shown that the proposed transitional breakage probability can 
be approximated by a series of ’power law’ models Kv^ with 
different exponents for different drop volume regions and 
that the value of the exponent increases as drop volume 
decreases, A functional relation for the rate constant K in 
the 'power law’ approximation of the transitional breakage 
probability is derived in terms of the parsimeters and physical 
properties of the system, A universal non-dimensionalized 
equilibrium drop volume distribution for agito/bed lean liquid- 
liquid dispersions is derived by analytical solution of the 
population balance equation simplified by order of magnitade 
estimates under the assumptions of constant coalescence 
frequency and uniform binary breakage. 

The prediction of equilibrium drop size distribution 
by solution of population balance equation necessitates the 
knowledge of breakage and coalescence. Even for lean 
dispersions, where the interaction between droplets is 
negligible, it is necessary to take coalescence into account 
which makes the population balance equation non-linear. Hence 
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it becomes necessary to make simplifying assumptions for the 
population balance eqLuation to be amenable to solution. 

Prediction of equilibrium drop size distribution from a 
statistical thermod 3 mamic -v-iewpoint overcomes this drawback as 
it considers only the total energy of the system i.e. the 
kinetic, potential and interfacial energies of droplets. A 
statistical thermodynamic model for equilibrium di'op size 
distribution is proposed. As the initeraction between droplets 
is negligible for a lean liquid-liquid dispersion, the system 
is assimied to consist of independent droplets, A canonical 
ensemble of droplets is considered for the derivation of 
equilibrium drop size distribution. The energy of a droplet 
is split into; (i) Kinetic energy of pulsating convective 
Liotion (ii) Kinetic and potential energies of oscillation and 
(iii) Potential energy of a droplet at its spherical equilibrium 
shape. A partition function for a droplet is derived, A 
droplet is found to have 6 degrees of freedom, A universal 
non-dim en si onalized drop volume density in terms of dimensionless 
drop diameter d/3!^2 derived from the partition function of 
a droplet. 

Prop volume distributions in a batch agitated lean 
liquid-liquid dispersion eventually approach an asymptotic 
distribution in the similarity variable (14-Kt)w^ if (i) the 
transitional breakage probability f (v) follows a ’power law’ 

Kv^ and (ii) the breakage is ’similar’ i.e, "tih® daughter 
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droplet size distribution depends only on the ratio of the 
daughter and parent drop volumes. The existence of such a 
similarity transformation is made use of to verify the ’power 
law' approximation for transitional breakage probability and 
also to obtain the exponent n, rate constant K and the 
daughter droplet distribution from the transients of drop 
volume distribution in a bcttch agitated lean liq.uid— liq.uid 
dispersion. 

Drop vol-ume distribution measurements were made in a 
jacketted glass vessel agitated by a 6 blade paddle impeller. 
Measurements were made at the dispersed phase volume fraction 
of 0,1 per cent in order to minimise drop coalescence.^ 

Micro encapsulation techniq.ue was employed to freeze the 
system at the time of measurement. This technique encap- 
sulates the droplets with a thin nylon film formed as a 
result of an interfacial polycondensation reaction. The volume 
distribution of drops was measured by means of electronic 
size analysis (Ooulter Counter), Special probe assemblies 
consisting of platinum electrodes with apertures of 2 mm 
and 0,8 mm diameter were used to monitor the droplets. The 
pulse signals from the probe were processed and fed to a 
Multichannel analyser to give numbei' count of droplets for 
different droplet volumes. This yielded volume distribution 
of droplets. 
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Transients of drop volme distribution in agitated lean 
liquid-liquid dispersions were measured for water - CCl^ + 
i-Octane (50 - 50 per cent) system at three different agitator 
speeds of 480, 420 and 300 vfith 3 inch diameter paddle 
impeller- and at the agitator speed of 480 RPM with 2 inch 
diameter paddle impeller. Transients of drop vol'ume distri- 
bution were also measured at the agitator speed of 300 RPM 
with 3 inch diameter paddle impeller for water-Anisole + CCl^ 
(80-20 per cent) and water-chlorobenzene systems. Experimental 
results showed that the transitional breakage probability 
could be approximated by a 'power law’ of varying exponents. 

The value of the exponents obtained from the experimental data 
agreed fairly well with the values obtained from the approxi- 
mation of the predicted transitional breakage probability by 
a series of 'power law' models of varying exponents. The 
computations of the transitional breakage probability as 
predicted by the model yielded verj^ large values of maximum 
stable drop volumes thus necessitating the introduction of a 
factor in the predicted transitional breakage probability in 
order to match its variation ’.-rith experimentally observed 
variation. It is believed that the introduction of this factor 
is necessary to accoxmt for Important inhomogeneities in the 
flow field of an agitated vessel. As the transitional breakage 
probability r(v) could be a,pproximated by a series of 'power 
law' models with different value of exponent for different 
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drop volinne range and for each of these regions Ev^*t is 
the similarity transformation at sufficiently large times, 
one "would expect that a transformation r(v)‘t would collapse 
all the drop volume distributions into a single curve 
irrespective of their region if the model for transitional 
breakage probability is indeed valid. The transformation 
erfo [ (2(21/3-1)^/362/3 2 , l/2(.)l/9 ^-5/18/j,o -1/3].^ 

seems to collapse all the drop volume distributions into a 
single curve thus vindicating the model for transitional 
breakage probability. 

The modified similarity transformation (^) -t^^^ 

V 

collapses all the drop volume densities in a region into a 
single curve thus vindicating "idle assumption of 'similar’ 
breakage made in the derivation of the similarity transformation. 
Hence the moments of daughter droplet distribution were 
computed from the moments of the modified drop volirme density 
f { (“) t^/^} and the daughter droplet distribution was 

' V 

estimated through orthonoimal polynomial expansion. The 
estimated daughter droplet distribution indicated that (i) the 
breakage mechanism seems to be the same for all systems and 
experimental conditions (ii) the normalised daughter droplet 
density is more or less same for all the cases and (iii) the 
daughter droplet density is broad. The functional dependence 
of the estimated rate constant K with respect to the agitator 
speed and interfacial tension agrees well with the model. 



The non-dim ensionalized equilibrium drop volume density 
in terms of dimensionless drop diameter (d/3^2) Predicted by 
the kinetic and statistical thermodynamic models agree very 
well with each other. It is very interesting to note that 
the equilibrium drop volume density arrived at throtigh two 
entirely different approaches predict more or less the same 
density. The experimental drop volume density fy(d/d^2) 
agreed fairly well with the gaussian approximation. However, 
there is a shift in the experimental drop volume density 
towards smaller dimensionless drop diameter (d/ff^^) when 
compared with the kinetic and statistical thermodynamic models 
and the agreement with the coalescence-redispersion model 
seems to be betteri Comparison of the three models with the 
experimental data of Chen and Middleman, Brown and Pit and 
Shinnar showed that the experimental data agreed fairly well 
with kinetic and statistical thermodynamic models and that the 
coalescenco-redispersion model did not agree with the experi- 
mental data. The discrepancy between the present experimental 
data of equilibrium drop volume density sind the experimental 
data of Chen and Middleman, Brown and Pit and Shinnar could' 
be due to loss of count of very small voltage pulses corres- 
ponding to droplets of very small volumes because of poor 
signal to noise ratio for these droplets. 
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®'inn “ Constant defined in equation (4.2.8) 

A = Work function 

"b = Constant defined in equation (4.2,8) 

mn 

'^mn “ Constant defined in equation (4.2.23) 

d = Drop diameter 

d ^2 “ Sauter mean drop diameter 

d = Drop diameter at which the equilibrium drop volume 

density exhibits a maximum. 

df = Area of a surface element 

df = Area of a spherical surface element 

s 

^mn ~ Constant defined in equation (4.2.23) 
d = Maximum stable drop diameter 

E = Internal energy 

E = Ensemble average energy of a droplet 

Emin = Minimum increase in specific interfacial energy of 
a droplet for fragmentation 

f = Eon-dimensionalized drop volume density in terms 

of droplet diameter 

f = Cimaulative distribution function in terms of 

similarity variable z 

fn = Normalised f 

f^ = Non-dimensionalized drop volume density in terms of 
droplet volume 

f ' = Drop volume density function in terms of similarity 

variable z 

F = Cumulative drop volume distribution function 

g = Cumulative volume fraction of daughter droplets 

expressed in terms of the ratio of daughter and 
parent drop volumes 
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g’ = Volume density of dau^ter droplets expressed in 
terms of the ratio of dau^ter and parent drop 
volumes. 

G = Cumulative distribution function of daughter droplet 

volume. 

H = Liquid depth in agitated vessel 

H = Hamiltonian of a droplet 

k = Boltzmann’s constant 

li = Rate constant in 'power law' expression for transitional 

breakage probability 

L = Impeller diameter 

M = Mass of a droplet 

n = Exponent in 'power law' expression for transitional 

breakage probability 

H = Speed of agitation, rotations/unit time 

= Ensemble average of total number of droplets 
= Weber number 

we p 

P = Density function of dau^ter droplet volume arising 

from breakage of a parent droplet in terms of number 

p = Momentum of a droplet 

p = Induced pressure distribution on the surface of a 

droplet 

p^ = Capillary pressure 

P = Power input 

P = Pressure 

(P.E) = Potential energy of a droplet at its spherical 

^ equilibrium shape 

^m ~ Associated Legendre function 

Pjj^ = Legendre polynomial of order m 
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Q = Partition function of a droplet 

r = Droplet radius 

t = Time 

T = Diameter of the mixing vessel 

T = Kinetic energy/- of oscillation 

T = Temperature 

T^ = Dispersion temperature 

u = Relative velocity between two points 

u^ = Characteristic velocity of an eddy of scale a 

V = Drop volume 

V = Maxomum stable drop volume 

V* = Drop volume at which the equilibrium drop volume 
density exhibits a maximum 

V = Dispersed phase volume fraction 

V = Volirme of liquid-liquid dispersion 
= Velocity of oscillation 

w = Frequency of fluctuation 

\Lm ~ natural frequency of oscillation of a droplet 
X = Hon-dimensionalized drop volume 

X = Variable representing the ratio of daughter droplet 

volume to parent drop volume 

X = Droplet diameter 

* 

X = Parameter defined in equation (4.3.2) 

2 = Similarity variable (1 + Kt)V^ 
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Greek Symbols ; 

a = Parameter defined in equation (3.2,12) 

®mn ~ Constant defined in equation (4.2.2) 

P = Moment of modified daughter droplet distrihution 

function for ’similar’ hreakage 

p = Parameter defined in equation (4.2,23) 

^mn “ Constant defined in equation (4.2.2) 

Y = Universal transitional "breakage probability in 

terms of non-dimensionalized droplet volume or 
diameter 

Y = Moment of daugh'te^ droplet distribution function 

for ’similar' breakage 

r = Transitional breakage probability 

X = Average number of eddies arriving on the surface 
of a droplet/unit time 

X = Moment of drop volume density in terms of modified 
similarity variable (v/v)t^/^ 

a = Interfacial tension 

g- = Standard deviation of equilibrium drop volume density 

p = Density of continuous phase 

e = Average power input/unit mass 

e = Energy dissipation rate/unit mass 

\i = Coalescence frequency 

p, = Moment of drop volume density in terms of similarity 

variable 

p = Viscosity of the continuous phase 

p, = Chemical potential 

= Mean number of droplets from breakage of a parent 
droplet 
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0 = Density function of droplet Tolniae in terms of number 

0 = Dispersed phase volume fraction 

f' = Velocity potential of the flow field inside a droplet 

= Modified velocity potential 
p = Generalized deformation coordinate 

5 = Non-dimensionalized droplet diameter 

^ = Deformation of a droplet from its spherical shape 

Subscr ipts 

X = X component 

y = y component 

2 = z component 

p = Moment-urn conjugate to p 

X = Momentum conjugate to X 



CHAPTER 1 


IHTROPUCTIOH 

Dispersed systems are of commoia occurrence in many 
industrially important processes. Dispersed phase reaction, 
liquid -liquid extraction, gas absorption and emulsion poly- 
merization are some such frequently encountered processes. 
Dispersed systems involve exchange of mass, momentum, energy 
or their combination between the dispersed and the continuous 
phases. The primary reason for dispersing one phase into the 
other in the form of droplets or bubbles is to increase the 
overall extent of such exchange through the increase of 
interfacial area. 

Agitated liquid -liquid dispersions are encountered 
mainly in stagevrise liquid -liquid extraction and emulsion 
polymerization. The interfacial area in such systems depends 
upon the drop size distribution, which is greatly influenced 
by drop breakage and coalescence. Moreover, the overall 
performance of an extraction stage or an emulsion polymeri- 
zation reactor is influenced by such factors as the drop 
breakage, the repeated coalescence and redispersion of drops 
in turbulent flow field, local mass transfer and reaction 
rates between drops and the continuous phase. Although a 
cohesive model of the whole cannot yet be formulated, it is 
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necessary for a designer to have some information on each 
of the factors mentioned above, in order to make an intelli- 
gent selection of design or scale-up rules for his apparatus „ 
The general approach to the calculation of rate 
processes in liquid -liquid dispersions has been through the 
definition of an overall transfer coefficient or reaction 
rate based on the total interfacial area and average properties 
(concentration and temperature). This approach ignores 
altogether the fact that the dispersed phase exists as a, 
population of droplets of different sizes and properties. 
Therefore, it does not take into account two important 
features of a dispersed system. Firstly, the transfer 
coefficient from an individual droplet directly depends on 
the size of the droplet which is influenced by breakage and 
coalescence; secondly, when the droplets coalesce, the 
gradients of concentration and temperature existing in the 
droplet before coalescence are destroyed and the resultant 
droplet attains the average concentration and temperature [1] , 
Moreover, the rate of breakage of 'primary droplets' [2] 
(Droplets formed due' to breakage of large globules) will be 
quite small. The droplet size range of these primary droplets 
would possess large surface to vol'ume ratio causing concen- 
tration and temperature changes at a rate likely to comp" 
with the breakage rate of these drops. Thus, it is necessary 
to take into account the simultaneity of droplet phenomena 
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and rate processes for prediction of overall rate and 
extent of transfer in liq.uid -liquid dispersions. 

The prohlem of prediction of transfer rates between 
dispersed and continuous phases in liquid -liquid dispersions, 
is therefore quite complex and must consider many factors 
that are not yet very well understood, Among these are the 
nature of drop breakage in a turbulent flow field with the 
associated formation of a wide distribution of drop sizes, 
the interaction between mass transfer and hydrodynamic effects, 
the repeated coalescence and redispersion of droplets in a 
dynamic environment. Though it is not possible to study all' 
these factors at a time due to their inherent complexity, 
it is believed that study of these factors separately will 
not only lead to a better ■understanding of the process but 
also will contribute greatly to the evolution of a rational 
methodology for the prediction of transfer rates in liq'uid- 
liquid dispersions. With this in view, the present study 
restricts itself to drop breakage in agitated liquid-liquid 
dispersions. Moreover, in case of lean liquid -liquid 
dispersion (low dispersed phase fraction) only breakage need 
be considered as coalescence will be negligible due to 
negligible interaction between droplets. 

Drop breakage in a turbulent flow field is random and 
hence is described by related probabilities, namely (i) the 
transitional breakage probability and (ii) daughter droplet 
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distribution [5]* -to attempt is made to relate these pro- 
babilities to the characteristics of turbulent flow field 
and physical properties of the system, A model is proposed 
for the transitional breakage probability of droplets exposed 
to a turbulent flow field based on the mechanism of droplet 
breakage due to their oscillation resulting from relative 
velocity and pressure fluctuations of the turbulent flow 
field. Using this model, population balance equation' was 
solved to predict the equilibrium drop size distribution. 
Transients of drop volume distribution in a batch agitated 
lean liquid -liquid dispersion was measured using encapsulation 
and electronic size analysis. Ramkrishna [4] proposed a 
method of evaluating the transitional breakage probability 
and daughter droplet distribution from the transients of drop 
volume distribution in a batch agitated lean liquid-liquid 
dispersion. This method was employed to calculate these 
probabilities from the experimental data which are then 
compared with the model, A model for equilibrium drop size 
distribution is also proposed from a statistical thennod3mamic 
viewpoint. This model considers the total energy of liquid- 
liquid dispersion but does not take into account explicitly 
drop breakage and coalescence. The measured equilibrium drop 
size distribution was also compared with the model. 

Chapter 2 gives a brief review of literature on 
drop breakage. Chapter 5 discusses the model for transitional 
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breakage probability and the evaluation of the probability 
and daughter droplet distribution from transients of drop volume 
distribution in a batch agitated lean liquid-liquid dispersion. 
Chapter 4 deals with the model on equilibrium drop size 
distribution from a statistical thermodynamic viewpoint. 

Chapter 5 discusses the experimental technique, apparatus 
and procedure. Experimental results and comparison between 
experimental data and the above models are presented in 
Chapter 6, Chapter 7 gives conclusion and recommendations. 



CHAPTER 2 


LITERATURE REVIEW 

Drop breakage in agitated liquid -liquid dispersions 
is complex due to the complex nature of the turbulent hydro- 
dynamic flow field to which the droplets are exposed. Since 
it is necessary to understand the droplet phenomena of breakage 
and coalescence in order to estimate the rate processes in 
agitated liquid -liquid dispersions, considerable attempts 
have been made to study drop breakage in agitated liquid- 
liquid dispersions, A brief review of the work accomplished 
in this direction is presented below. 

2.1 DROP BREAKAGE ME CHARI SMS; 

In liquid-liquid dispersions most of the ordinary 
dimunitive processes function either by injection of one phase 
into the other or by inducement of turbulence in the continuous 
phase. In case of injection, the potential energy of the 
phase to be dispersed is converted into kinetic energy, and at 
the same time usually interaction takes place with the conti- 
nuous phase. This interaction generates forces in the phase 
to be dispersed which may result in its breaking up. In case 
of turbulence, it is the kinetic energy of the turbulent 
motion in the continuous phase that brings about the breakup 
of the other phase. 
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In general, the disintegration process takes place 
in stages. When one fluid is just squirted into another or 
•when agitation has just started, the phase to he dispersed 
is at first present in hulk. The deformation of this hiilk 
of fluid and its initial breakup into chunks of fluid, which 
may breakup further into smaller parts^has been much studied 
[5,6], The mechanism of breakup of the bulk of fluid into 
globules seems to be the penetration of lamellae and ligaments 
of one fluid into another [6], These ligaments then breakup 
into globules, which may further split into smaller parts. 
Important contributions have been made by theoretical and 
experimental studies of the instability of liquid films and 
cylindrical liquid ligaments [6,7], These should give insight 
into the rupture of the lamellae into glob-ules. These theories 
can be applied to some simple types of atomisation. But, in 
most cases, the breakup of bulk of the dispersed phase into 
globules is only the first stage of disintegration. This is 
usually followed by further breakup of these globules into 
smaller drops. When the viscosity and interfacial tension of 
the dispersed phase are not very high, this first stage of 
disintegration, viz. disintegration of bulk of liquid into 
globules, is extremely rapid and chaotic [8], Since this 
first stage of disintegration is extremely rapid, it generally 
does not eontrol the size distribution of droplets formed 
after further disintegration of globules. 
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The second stage of breakage is the further 

disintegration of the globules into smaller droplets as 

stated earlier. The study of this stage of disintegration 

is very important as it seems that this 'Will control the 

sizes of droplets eventually formed in a breakage process. 

Fundamental work on the splitting of drops in viscous flow 

was done by Taylor [9] and Tomotika [10], Taylor [9] 

developed a theory based on a few obvious assimaptions , such 

as, small deformations, no slip at the interface and continuous 

shear stress across the interface. He considered two types 

of flow, namely, the couette flow and the plane hyperbolic 

flow. According to Taylor’s theory^ the deformation of drops 

is predominantly determined by the generalised Weber group 
\i S D 

%e “ ”"^0 * where S is the maximum velocity gradient in the 

flow field of the continuous phase. He determined the value 
of S at which breakup of a drop occurs, that is 
He observed that the droplet exposed to the shear field 
elongated into a threadlike form. The disintegration of this 
thread led to formation of drops which are of the order of 
l/lOOth of the size of the original drop. Ramscheidt and 
Mason [11] studied the deformation and burst of fluid drops 
in a liquid subjected to both hyperbolic and shear flow for 
a wide range of viscosity ratio. At low deformation, excellent 
agreement with the theoretical equations of Taylor was found. 
The critical conditions of burst were also found to be in 
reasonable agreement with the theory. 
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Hinze [8] studied the brealcup of globules in viscous 
flow, brealcup of drops exposed to airstream and emulsification 
in turbixLent flow. He proposed that the deformation of a 
globiiLe can be due to viscous stress or dynamic pressure and 
that the interfacial tension will give rise to a surface 
force that will in general counteract the deformation. He 
identified two dimensionless groups which would control the 
deformation and breakup: a Weber group and a viscosity 
group ^Yi* proposed that breakup occurs when 

exceeds a critical value • Experimental observations 

showed that the mechanism of breakup was found to be quite 
simple when is equal to or sli^tly higher than Qpj_.^ • 

The more exceed? its critical value, the more complicated 

this mechanian becomes. For Qrit mechanism 

is very complex and the disintegration process is more or 
less chaotic [8] . This is probably the reason why the second 
stage of breakage is simple and amenable to study compared 
to the first stage of disintegration, Hinze [8] also showed 
that in case of droplets exp)osed to an airstream and emulsi- 
fication in turbulent flow, dynamic pressure forces rather 
than viscous forces were responsible for deformation and 
breakup. In case of emulsification in turbulent flow, the 
dynamic pressure forces were caused by changes in velocity 
at most equal to the diameter of the drop [8], 
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2.2 HTDRODYFMIC CHARACTERISTICS OR THE RLO¥ RIELD IN AN 
A&ITATED VESSEL: 

In order to study the drop "breakage in agitated 
liq.uid -liquid dispersions, it is necessary to have knowledge 
about the hydrodynamic characteristics of the flow field 
in an agitated vessel since the interacting forces responsible 
for defoliation and eventual breakup of a droplet will depend 
on the nat'ure of the flow field. Considerable amount of 
work has been done regarding the flow pattern, power consumption 
and the use of dimensional analysis in agitated vessels. 
Excellent reviews of the voluminous literature on flow fields 
in agitated vessels have been made [12, 13]. 

The flow regions of a stirred tank may be classified 
into impeller stream flowing radially towards the wall, the 
recirculation regions and the region surrounding the impeller 
where the recirculation loops meet. An extensive study of 
the flow pattern in a 60 cm tank has been reported by Nagata 
et al. [14]. They divided the recirculation flow pattern 
into roughly four zones; (i) the horizontal discharge jet, 
running from the impeller to the tank wall (ii) the vertical 
flow along the wall, (iii) the horizontal flow returning the 
fluid to the vicinity of the tank axis and (iv) the vertical 
flow along the axis towards the impeller. Power consumption 
in stirred vessels has been measured extensively and reported 
in the form of correlation for dimensionless power number as 
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a fimction of Reynolds numler for different types of 
impellers [15, 16], Measurement of the discharge rate of 
impellers has also "been reported as a function of agitator 
speed and diameter [17], 

Considerable attempts have recently been made to 
understand the nature of turbulence in agitated vessels. 

Cutter [18] studied turbulence in a batch baffled stirred 
tank with a disc turbine by means of a photographic techniq.ue. 
Most of the results presented were about the impeller stream. 
The important turbulence parameters obtained were intensities 
of turbulence, macroscales in the radial and tangential 
directions and estimates of energy dissipation rates in various 
parts of the tank from the energy equation. This work 
indicated that half of the total energy dissipation in the 
tank occurred in the impeller stream, a fifth in the impeller 
region itself and the remaining in the rest of the tank 
exhibiting important inhomogeneities in the flow regime. An 
important conjecture of Cutter [18] was the presence of 
intermittancy in the impeller stream, Schwartzberg and 
Treybal [19] used a photographic technique to determine the 
mean and fluctuating velocities in a baffled agitated vessel 
with a disc turbine. The measurements were made outside the 
impeller stream, Rim and Manning [20] presented radial 
intensity spectra as measured with a transducer probe. A 
region of slope -5/5 was observed in their spectra. 
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Muzumdar et al, [21] measured the turhulence parameters in a 
turbine stirred baffled tank with air as the medium using 
hot wire anemometer in the impeller stream. They observed 
large periodic components in the velocity near the impeller. 

They could not detect the region of -5/3 slope in the energy 
spectrum due to superimposition of periodic non-random velo- 
cities on the random velocity fluctuations. They also measured 
the probability density of velocity fluctuation which was 
found to be far from gaussian. I. Komasawa, R. Kuboi et al, [22] 
studied the turbulence in turbine stirred baffled vessel 
using tracer particles with the help of photography. Energy 
spectrum obtained by fourier analysis of autocorrelation 
showed a region of -5/3 slope, M. Ananda Rao and Broadkey [23] 
measured the mean velocities, EMS value of velocity fluctu- 
ation, one-dimensional energy spectra, normalised autocorre- 
lation and probability density of velocity fluctuation in a 
continuous turbine stirred vessel with water as the medium 
using hot wire anemometer. They also observed periodic 
velocity components near the impeller. One dimensional energy 
spectrum exhibited a region of —5/3 slope over a decade of 
wave number range. The probability density of velocity 
fluctuation was found to be skewed negatively and far from 
gaussian due to intermittancy and non-isotropic conditions. 

A, A, G-tlnkel and M.E, Weber [24] measured the turbulence 
parameters in a baffled stirred tank agitated by a six blade 
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disk turbine with air as the medi-um both in the impeller 
stream as well as in the bulk of the tank using hot wire 
probe. They also observed periodic velocity components in 
the impeller stream near the impeller. Their one dimensional 
energy spectrim in the impeller stream consisted of three regions, 
namely, the region of zero slope followed by the regions of 
-5/5 and -7 slopes. The nature of the energy spectrum was 
very much similar to the spectrum obtained by others [20, 21, 

22, 23]. They calculated the rate of energy dissipation 
within the impeller stream by using the energy balance equation. 
The calculations showed that very little energy is dissipated 
within the impeller stream indicating that most of the energy 
dissipation took place in the bulk of the tank. This was 
in contradiction with the findings of Cutter [18], A, A, GtLnkel 
and M,E. Weber [24] claimed that the inhomo geni ties observed 
by Cutter [18] were due to experimental errors in the measure- 
ment of velocity fluctuation. The measurements of A. A, GtLnkel 
and M,E, Weber [24] in the bulk of the tank showed that the 
turbulence in the biilk of the tank is isotropic and the 
velocity fluctuations did not change appreciably with position. 
One - diieensional energy spectrum obtained in the bulk of 
the tank also showed regions of zero, -5/5 and -7 slopes. The 
calculations of energy dissipation rates from the one- 
dimensional energy spectrum seem to confirm the view that 
most of the energy was dissipated in the bulk of the tank. 
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They compared their value of macroscale with the values 
obtained by Cutter [18], Rao and Broadkey [-23] and Chao et al. 
X25] and showed that the turbulence in the bulk of the tank in 
the -5/3 slope region of the energy spectrum was independent 
of the viscosity of the medium and depends only on s, thus 
vindicating Kolmogorov's theory of local isotropy, 

2,3 experimektai iftestigatiors or drop breakage in agitated 

LIQUID-LIQUID DISPERSIONS: ' 

The study of agitated liquid -liquid dispersions has 
been the subject of numerous investigations [26, 27]. In such 
systems, the drops of the dispersed phase are continuously 
broken and coalesced at a rate depending on the local hydro- 
d 5 mamic characteristics of the turbulent flow field and the 
physical properties of the system. Earlier experimental 
investigations on agitated liquid -liquid dispersions focussed 
on measurements of average droplet diameter by the use of 
light transmission, photographic techniques or chemical methods. 
Such experimental measurements resulted in correlations for 
average drop diameter in agitated liquid-liquid dispersions. 

Drop size distribution measurements were also made with 
photography, encapsulation, sedimentation and electronic size 
analysis, A niambdr of investigators reported experimental 
data on interfacial area and drop size distribution in flow 
systems. Drop .breakage in agitated liquid-liquid dispersions 
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has been analysed by many using Kolmogorov* s theory of local 

isotropy. The applicability of the above theory in case of 

agitated liquid -liquid dispersions was illustrated with 

experimental evidence by Shlnnar and Church [ 28] , They 

characterized breakage in turbulent flow field by a maximum 

stable drop diameter and derived an equation for the same as 

a function of impeller diameter, agitator speed, density and 

interfacial tension. Similar equations were derived by Hinze[8], 

levich [29] and Chen and Middleman [30], For drops much larger thai 

the microscale of turbulence t], consideration of the balance 

between inertial and interfacial forces leads to a maximum 

stable drop diameter d This has been confirmed 

max 

experimentally for a variety of dispersed phase systems, covering 

a wide range of physical properties and experimental techniques 

[28, 30, 31, 32, 33, 34]. For drops smaller than t}, viscous 

effects predominate, and Taylor’s analysis [9] predicts 

d^ax“ *^[33] • Experimental evidence of this dependence 

has been reported for gas-liquid systems [31], and for low 

interfacial tension systems [34], although in both cases the 

drop size was greater than the microscale of turbulence. 

Vermuellen et al. [31] measured the interfacial area 

and mean drop diameter for agitated liquid -liquid dispersions 

using light transmission technique and arrived at a correlation 

—0 6 

for mean drop diameter which showed a dependence of * • 

The form of this correlation was found to be identical to the 
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theoretical equation derived by Shinnar and Church [28]. 

However, the correlation showed a dependence of mean drop 
diameter on dispersed phase fraction 0 which indicated that 
coalescence was also important in their experiments. They 
also studied the effect of and A p on mean drop diameter 

and found that they were relatively unimportant. Their Weber 
number was based on the mean density of the dispersion. 

Their experimental data exhibited deviation of as much as 
+ 60 per cent from the correlation. Only limited geometries, 
two tanks of different sizes and one impeller were employed 
in their experiments. Calderbank [32] measured the inter- 
facial area in agitated liquid -liquid dispersions with concen- 
tration of the dispersed phase ranging from 0-20 per cent. 

The variation in the physical properties of the systems studied 
was limited. Their correlation for the sauter mean drop 
diameter d ^2 showed a dependence of Ihe maximum 

deviation was found to be + 10 per cent. Rodger et al. [35] 
measured average drop size in agitated liquid -liq'uid dispersions 
using equal volumes of both liquids. Their correlation for 

showed a dependence of This agreed quite well 

with the relation for minimum stable drop diameter for coales- 
cence of drops derived by Shinnar [35] > using the concepts of 
local isotropy. This is understandable, as in Rodger et al.’s 
[35] experiments coalescence may be dominant due to high 
dispersed phase fraction. Chen and Middleman [30] measured 
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the equilibrium drop size distribution for agitated lean 
liquid-liquid dispersions. Their experimental sauter mean 
drop diameter ^^2 agreed well with the correlation of Yermuellen 
et al . [313. Chen and Middleman [30] established the normality 
of the equilibrium drop volume distribution and showed that 
the equilibrium drop volume distributions for different 
experimental conditions and systems can be represented by a 
single distribution when reduced in terms of dimensionless 

•i 

drop diameter d/J^2« Shinnar [ 33 ] measured drop size distri- 
butions of molten wax dispersed in hot water using photo- 
micrography* He showed that it is possible to obtain turbulence 
stabilized dispersions under certain conditions if the maximum 
stable drop diameter for drop breakage is above the minimum 
stable drop diameter for coalescence. His experimental data 
on equilibrium drop volume distributions when reduced to 
dimensionless drop diameter d/J^2 seem to agree fairly well 
with the correlation of Chen and Middleman [50], Brown and 
Pit [36] measured equilibrium drop size distribution in kerosene- 
water system using photography for low dispersed phase fractions. 
The experimental data suggested that coalescence was absent. 

The equilibrium drop volume density in terms of dimensionless 
drop diameter d/d^2 agreed fairly well with the correlation 
of Chen and Middleman [30], Observation of an equilibrium 
drop size distribution in batch lean liquid -liquid dispersions 
[30, 33 f 36] are indicative of restricted form of breakage 
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that lends credibility to the concept of maximtan stable drop 
diameter. Since coalescence rates are expected to be negligible 
for a lean liquid -liquid dispersion due to negligible inter- 
action between droplets, equilibrium drop size distribution 
in a batch agitated lean liquid-liquid dispersion does in fact 
correspond to a sluggish dynamic equilibrium between breakage 
and coalescence with the equilibrium density centred more or 
less around the maximum stable drop diameter [2]. Sprow [ 34 ] 
measured drop size distributions in the vicinity of a turbine 
impeller for lean iso-octanc dispei’sion in salt water using 
electronic size analysis (coulter counter). The experimental 
results of d agreed well with the correlation of Vermuellen 
et alo [31]. Sprow [37] also measured drop size distributions 
in a strongly coalescing agitated liquid -liquid system. The 
drop size was found to vary considerably with position and 
drop sizes were found to be much smaller near the impeller ■ 
and bigger outsrde the impeller stream. Such strong variation 
of drop size with position was also observed by Vanderveen [38] 
and Cald erb ank [32]. The experimental results of Sprow [37] 
indicated that drop breakup predominated near the impeller, 
whereas coalescence of drops was controlling in other locations. 
However, the experimental measurement of drop sizes in water- 
ed^ + i-octane system at low dispersed phase fractions using 
micro encapsulation by Mylnek and Resnick [39] showed that 
there is almost no dependence of drop size on location. This 
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was due to the fact that coalescence rate was low comp3,red 
to the circulation time. Comparison of mean drop diameters 
obtained using different impellers showed that the criterion 
of eq.ual power per unit ■volume cculd be used for estimating 
drop size when going from one geometry to another, not too 
different, geometry at moderate impeller speeds. Park and 
Blair [ 40 ] studied the droplet interaction phenomena in lean 
agitated li<luid— liquid dispersions using cine photography and 
light transmission technique. They observed coalescence and 
breakup events and measured the rate of coalescence at various 
positions in the tank. They found that drop breakup occurred 
near the impeller and the droplet coalescence predominated at 
other locations and that beyond distances from the impeller 
region of order of only l/6th the impeller diameter, breakup 
was virtually non-existent, 

Yarious investigators have reported average drop size 
and drop size distribution measurements for continuous flow 
agitated liquid-liquid dispersions. In continuous flow systems, 
drop size will be determined by the residence time of the 
dispersed phase in addition to the hydrodynamic conditions 
and physical properties. No correlation has been developed to 
acco'unt for this. Schindler and Treybal [ 41 ] and Co'ulaloglou 
and Tavlarides [ 42 ] observed that drop size varied with 
position in a continuous flow system, CoiiLaloglou and 
Tavlarides [ 42 ] measured the drop size distributions in the 



20 


circulation and impeller regions. They also measured coales- 
cence frequencies. Differences observed in the size distri- 
butions between the regions were found to be smaller than 
would be expected considering the largely different hydrodynamic 
characteristics of the flow field. They attributed the relative 
spatial homogeneity of the dispersion to the low coalescence 
frequency amongst the droplets relative to the average internal 
circulation frequency of the fluid elements from the impeller 
tip and back and suggested that the ratio of the coalescence 
frequency to the circulation frequency can serve as a measure 
to determine whether or not one shoiiLd treat a dispersion as 
being spatially homogeneous. It was observed that the dispersion 
was homogeneous as this ratio approached zero. 

2,4 MODELLING OF DROP BEBAIIAGB IN AGITATED LIQUID-LIQUID 

DISPERSIONS J 

The general mathematical formulation for the description 
of particulate systems was put forward by Hulbert and Katz [43] 
in which they describe the particulate system by a set of phase 
coordinates viz, the spatial coordinates of the particle and 
the properties of the particlOi They developed a partial 
differential equation for the number density of particles in 
the phase space. Curl [44] employed population balance 
approach to describe dispersed phase mixing in agitated liquid- 
liquid dispersions, A simplified model was employed in which 
only uniform drop sizes were considered. The model assumed 
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two drop coalescences, the prohability of coalescence being 
the same for every droplet. According to the model, the 
coalesced drops redisperse immediately to foim new drops of 
equal size and coalescence and redispersion phenomena have no 
effect on mass transfer. The model is unrealistic as it does 
not include effects due to a distribution of drop sizes, 
uneven breakup of drops or a distribution of probabilities of 
coalescence with drop size or concentration. Baynes and 
Laurence [45] employed a trivariate distribution of droplet 
mass, amount of solute and age and population balance equation 
to describe mass transfer in agitated liquid -liquid dispersions. 
They considered homogeneous batch stirred tanks with equal 
rates of coalescence and breakage and with only binary coales- 
cences followed by instantaneous breakup to form two drops of ' 
equal volume. They considered the simultaneity of droplet 
phenomena and mass transfer by introducing the concept of 'age' 
of a droplet and by employing a multivariate density for the 
calculation of mass transfer rates, Valentes, Bilous and 
Mundson [3] developed a comprehensive mathematical model for 
drop breakage in agitated liquid-liquid dispersion and predicted 
the drop size distribution evolved due to breakage by solving 
the population balance equation. They represented drop breakage 
by related probabilities, namely, the transitional breakage 
probability, daughter droplet size distribution and average 
number of daughter droplets formed upon breakage. They 
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considered limited and nnlimited "breakage of drops. They 
assumed a ’power law’ model for transitional breakage probability. 
The model for limited breakage included the effects of agitator 
speed and physical properties of the system by taking into 
account the maximum stable drop diameter. Their calculations 
showed that the drop size distribution was sensitive to the 
average breakage time and the average number of daughter 
droplets formed upon breakage. Fon-ideality in mixing vessels 
was also found to affect the drop size distribution in a 
continuous system. Valent es and imundson [46] predicted the 
drop size distribution for a coalescing agitated liquid-liCLuid 
dispersion by proposing a model for coalescence frequency. 

The calculation of transients of drop size distribution in a 
batch agitated liquid-liquid dispersion with limited breakage 
indicated the existence of a dynamic equilibrium between 
breakage and coalescence for large times with the equilibrium 
drop size density independent of initial size distribution, 

Bajpai, Ramkrishna and Prokop [47] proposed a coalescence- 
redispersion model for equilibrium drop size distribution for an 
agitated liquid-liquid dispersion. They assumed that drop 
coalescence is followed by immediate redispersion into two 
drops whose sizes are distributed according to a uniform 
distribution in the ratio of the daughter drop' volume to 
parent drop volme. They solved the pop^llation balance equation 
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under the assumption of constant coalescence frequency and 
obtained an exponential distribution for equilibrium drop 
■volume distribution. It was found that the prediction was 
good for flat blade t'urbine agitated systems, lor circulation 
stirring systems, the model predictions were substantially 
improved by assuming that the coalescence frequency is 
proportional to the diameters of the interacting drops. Steidl 
[48] proposed a statistical thermodynamic model for eq'uilibrium 
drop size distribution in agitated liquid -liquid dispersion. 

The equilibrium drop size distribution was expressed in terms 
of a new temperature scale. However, the coordinates and the 
corresponding momentum coordinates of the defined phase space 
for liquid -liquid dispersions were not conjugate to each 
other. Although not directly applicable to agitated liquid- 
liq'uid dispersions, the work of Collins and Knudson [49] is of 
interest. They measured photographically the drop size 
distributions produced by turb'ulent pipe flow of two immiscible 
fluids and developed a stochastic model to describe the 
breakage of drops. They assumed that coalescence is absent. 
They incorporated maximum stable drop diameter d^„_ in their 
stochastic model. They also introduced a parameter d^ ^ which 
referred to the drop size above which drop breakup occurred 
instantaneously. One of the many post'ulated stochastic models 
was found to represent the experimental data fairly well. 

The introduction of the parameter d^ ^ corroborates the view 
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that drop breakage occurs in two stages and that the first 
stage of breakage, Tiz. breakage of large lumps of fluid, 
is instantaneous, Shah and Ramkrishna [1] proposed a 
population balance model for mass transfer in lean liquid -liquid 
dispersions. They accounted for the simultaneity of drop 
breakage phenomena and mass transfer by using a trivariafe‘= 
density in terms of drop size, concentration and age. They, 
used a 'power law' model for transitional breakage probability 
similar to the one proposed by Valentes, Bilous and imundson[3]. 
They assumed breakage to be limited and introduced a parameter 
akin to d^ ^ introduced by Collins and Khudson [49], in their 
model for transitional breakage probability, Zeitlin and 
Tavlarides [50] developed a model which accounts for complex 
fluid-fluid interactions and hydrodynamic effects in an agitated 
liquid-liquid dispersion. The mixing vessel was divided into 
different regions to account for the hydrodynamic inhomogeneities. 
Based on the model, the system was simulated using Monte-Carlo 
simulation technique. Both micromixing and macromixing effects 
on particle size distribution were taken into account by 
simulating the particle movement in the vessel and interaction 
between the droplets. The dispersed phase was ass-umed to 
consist of spherical drops with non-uniform size distribution 
and only biparticle coalescences were assumed to occur. 

Drop breakage was assumed to be limited and binary. Maximum 
stable drop diameter and a parameter d^ ^^ analogous to 

d^ Q introduced by Collins and Khudson [49] were incorporated 
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in the breakage model. Probability of breakage of a droplet 
■was ass-umed to follow a ’power law’. The parameters of the 
model were estimated by comparing the calculated and experi- 
mental distributions. 

Swift and Freidlander [51] in their study of coagu- 
lation of hydrosols have introduced the concept of self- 
preserving distributions which depends upon the existence 
of a similarity transformation between particle size variable 
and time, Eapur [52] has provided an application of this 
concept to problems in size reduction. Pillipov [55] in 
an interesting study of batch particle splitting processes 
has shown the existence of explicit similarity transformation 
■under the conditions of ‘power law’ model for transitional 
breakage probability and 'similar' breakage, Ramkrishna [4] 
has shown the method of estimation of the 'power law' constants 
of transitional breakage probability and daughter droplet 
distribution from the transients of drop volume distribution 
in a batch agitated lean liquid-liquid dispersion, A 
preliminary analysis of the experimental data of transients 
of drop size distribution in batch agitated lean liquid-liquid 
dispersion obtained by Madden and McCoy [54] showed that 
'power law' model for transitional breakage probability is 
indeed valid, Co^ulaloglou and Tavlarides [55] proposed 
phenemenological models for drop breakage and coalescence 
in a turbulently agitated liquid-liquid dispersion, Pop'ulation 
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balance eq.uation with, the proposed models for breakage and 
coalescence was solved using iterative technique to predict 
the evolved drop size distribution. Drop breakage was 
described by a transitional breakage probability, A model for 
transitional breakage probability was proposed using Kolmogorov’s 
theory of local isotropy. It was assumed that the fraction 
of drops breaking due to dynamic velocity and pressure 
fluctuations is proportional to the fraction of turbulent 
eddies colliding with the droplets that have a turbulent 
kinetic energy greater than the interfacial energy pf the 
droplet. This was estimated under the assumption that the 
probability distribution of kinetic energy of eddies is 
exponential. Concept of breakage time was introduced in the 
definition of transitional breakage probability. Breakage 
time was estimated by assiming that the motion of the centres 
of mass of the daughter droplets to be formed is similar to 
the relative motion of two lumps of fluid in a turbulent 
flow field. The model was found to predict correctly the 
effect of impeller speed on maximum stable drop diameter. 



CHAPTER 3 


A MODEL FOR TRANSITIONAL BREAKAGE PROBABILITY OP DROPLETS 
IN AGITATED LEAN LIQNID-LIQUID DISPERSIONS 

Drop si25e distributions in liquid -liquid dispersions 
evolve due to the dynamics of breakage and coalescence. 
Evolution of drop size distribution is important for prediction 
of, interfacial area, heat and/or mass transfer in liquid- 
liquid dispersions. Attempts have been made [1,3,46] to 
predict the drop size distribution in liquid -liquid dispersions 
using population balance equation# If considerations are 
restricted to lean dispersions (i.e. with low dispersed phase 
fractions) only breakage need be considered as coalescence is 
negligible due to very little interaction of droplets. Drop 
breakage in an agitated liquid-liquid dispersion is random 
due to random velocity and pressure fluctuations responsible 
for deformation and breakup of drops# Dynamics of drop 
breakage in agitated liquid -liquid dispersions has been 
characterised by transitional breakage probability r(v) and 
daughter droplet size distribution G(v,v' ) [3,4] defined as, 

r(v)dt = Probability that a drop of volume v will 
break in time interval t and t+dt# 
and G(v,v*) = Volume fraction of daughter droplets with 

volume less than v formed from breakage of a 
droplet of volume v’ # 
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Quantitative prediction of droplet dynamics in lean 
dispersions requires information regarding transitional 
breakage probability r(v) and daughter droplet distribution 
G(v,v'), A model for transitional breakage probability of 
droplets in agitated lean liquid -liq-uid dispersions is 
proposed based on the mechanism of breakage of droplets due 
to their oscillations resiilting from relative velocity fluctu- 
ations. A universal non-dim ensionalized equilibrium drop 
volume distribution for batch agitated lean liquid-liquid 
dispersions is derived by analytical solution of population 
balance equation simplified by order of magnitude estimates. 
Since it is not possible to make a direct experimental measure- 
ment of transitional breakage probability P (v) , one has to 
devise an indirect method of estimation of r(v) in order to 
verify the model, Ramkrishna [4] has shown the method of 
estimating 'power law' constants of transitional breakage 
probability and daughter droplet size distribution from the 
transients of drop volume distribution in a batch agitated 
lean liquid-liquid dispersion. The details of the method are 
also presented. 

5.1 MODEL EOR TRANSITIONAL BREAKAGE PROBABILITT: 

A droplet e^qposed to a tuurbulent flow field will be 
subject to both inertial as well as viscous forces [55] • If 
the droplet is much larger than the microscale of turbulence, 
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the viscous forces caai be neglected. In agitated vessels, 
the Reynolds numbers are so high that the microscale of 
turbulence is found to be much smaller than the size range 
of droplets encountered [28], In such a case, the droplet 
will oscillate about its spherical equilibrium. The osci- 
llations are brought about by the kinetic energy of the 
turbulent motion in the continuous phase [28, 33], or 
equivalently, by the relative velocity fluctuations between 
points near the vicinity of the droplet surface. In other 
words, the kinetic energy of the turbulent motion brings 
about the increase in the interfaciaJ. energy of the droplet 
through deformations. Fragmentation of the droplet occurs, 
if the turbulent motion provides the minimum increase in 
the interfacial energy for breakup. 

In a turbulent flow field, velocity fluctuations at 
a point can be thought of as due to arrival of eddies of 
different scales (frequencies). Similarly, the relative 
velocity fluctuations around the surface of a droplet exposed 
to a turb-ulent flow field can be viewed as due to arrival 
of either an eddy or eddies of different scale (frequency) 
on the surface of the droplet. 

At this point, we compare the time scales of (l) arrival 
of an eddy on the surface of a droplet and (2) oscillation 
of the droplet. The time scale of oscillation of a droplet 
should correspond to the frequency of the eddy causing these 



30 


oscillations. Since oscilla,tions can be caused only by 
eddies of scale smaller than the diameter of the droplet, 
for sufficiently small droplets the time scales of oscillation 
should be very small. Moreover, oscillation of a droplet 
around its spherical equilibrium position necessitates the 
time scale of arrival of an eddy to be much larger than the 
time scale of oscillations, for otherwise, the effect of 
an eddy will not pei-sist long enough for the oscillations 
of a droplet around, its spherical equilibrium position to 
occur. If the time scales of arrival of an eddy is comparable 
to or less than the time scale of oscillation of the droplet, 
the effect of eddies will be to cause an irregular or bulgy 
deformation of the droplet, Hinze [8] considered bulgy 
deformation as one of the possible mechanisms of deformation of 
a globule exposed to a turbulent flow field. 

At sufficiently high Reynolds numbers, the flow field 
in an agitated lean liquid -liquid dispersion is locally 
isotropic [20, 21, 23, 24, 28], Hence we need to consider 
only the relative velocity fluctuations between any two 
diametrically opposite points on the surface of the droplet [33] 
i.e, the oscillating droplet can be considered as a one- 
dimensional simple harmonic oscillator. The points A and B 
represent any two diametrically opposite points on the droplet 
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surface. The restoring action of the drop due to inter- 
facial tension can he represented hy the spring connecting 
the points A and B, Since the size range of droplets is 
much above the microscale (unless the agitation speeds are 
very high), viscosity effects are negligible. Therefore, 
the simple harmonic oscillator can be considered undamped, 
Bor an undamped oscillator, the sum of potential and kinetic 
energy is constant, the potential energy signifying the 
interfacial energy of the droplet. Taking the datum of 
interfacial energy to be tha.t at the spherical equilibrium 
position. 


Maximum interfacial energy = Maxim-um kinetic energy of 

oscillation 


Thus, 


Maximum specific inter- 
facial energy 

Maximum increase in the 
specific interfacial 
energy above the 
equilibrium value 


= Energy of the arriving 
eddy 

= Energ7 of the arriving 
eddy 


Here we make a postulate that the number of eddies 
arriving on the surface of a droplet is a poisson process 

i.e. 

(i) The numbers of eddies arriving on the surface 
of a droplet in disjoint time intervals are 
independent. 

(ii) Pr [one eddy arrives on the surface in time 
interval t, t+dt] = X dt + O(dt^) 
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(iii) Pr [more than one eddy arrives on the surface 
in time interva,! t, t+dt] = 0 (dt^) 
where X = average number of eddies arriving on the surface 
of a droplut/unit time. 

Since we need to consider only the inertial subrange 
of Eolmogorov* s imiversal spectrum and the effect of eddies 
of scale less than or of the order of microscale of turbulence 
is negligible, we can neglect the effect of |i. on X . However, 
we can expect , the parameter of the Poisson process, to 
depend both on the diameter of the droplet and power input 
per unit mass e. In the absence of a priori information 
about the dependence of X on the droplet diameter and e, we 
a,ssume it to be independent of droplet diameter and e i.e. X 
can be regarded as a constant. 

We consider a small time interval r, and define two 
events A and B as 

A ; in eddy arrives on the surface of a droplet in 
the time interval v, 

B i Ihe arriving eddy has energy greater than or 

eq.ual to the minimum increase in the interfacial 
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Pr [Droplet of voItjuig v breaks in time interval t] 

= P(A) . P(B j A) 

= XT • P(B I A) 

Therefore, 

D'fc _ Prf Droplet of volmne v breaks in time interval t1 
% 0 

= r(v) = X • P(B j A) (5.1.1) 

The interfacial energy of a droplet of volume v 

= a ,.1/3 62/3 v2/3 

It can easily be shown that the increase in the interfacial 
energy re(iuired for fragmentation is minimum if binary eq.ual 
breakage occurs. 

Therefore, minimum increase in the interfacial energy for 
fragmentation 

= (2^/^ -1) cy tt; 1/3 g2/3 ^ 2/3 

This energy is provided by the kinetic energy of oscillation 
of the droplet. 

Therefore, minimum relative velocity of oscillation requiired 
for fragmentation is given by, 

I ( P v) = (2l/3_i) g g2/3 ^/3 (3.1.2) 

i 7 (3.1.3) 

If the characteristic velocity of an eddy of scale a is u^, 

o 

then its kinetic energy is ^ and for breakage to occur, we 

must havs | ^ i u^^^. 
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In order to evaluate r('v), we need to know the 
prohahility distribution of the relative velocity between 
two points separated by a distance r. Experimental measure- 
ments of probability distributioh of velocity fluctuation at 
a point in the impeller stream in an agitated vessel [21,23] 
have shown that this probability distribution may not be far 
from G-aussian, The gaussian approximation may be better at 
higher impeller speed (or higher Reynolds number) as inter— 
mittant flow conditions and non-isotropic conditions tend 
to become less pronounced, G-iven that the probability distri— 
bution of velocity at a point in an agitated vessel is 
normal, the most natural assumption about the joint probability 
distributions of velocities at two points is that it is 
normal. This, according to Batchelor [56] follows from the 
central limit theorem, since the turbulent velocities are 
subject to the influence of a large number of random eddies 
or local flow. On the other hand, such an assumption cannot 
be exactly true especially for very small values of r, 
because it amounts to neglecting non-linear terms in the 
dynamical equations, or, equivalently, the non-linear 
modulations of the energy spectrum [56], Due to lack of 
adequate information regarding the joint probability distri- 
bution of 2-point velocities in agitated vessels, we make 
the assumption that the joint probability distribution is 
normal. Prom this assumption, it follows that the probability 
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distribution of relative velocity between two points is also 
normal. 

If u is the relative velocity between two points 
separated by a distance d at any time t, from Kolmogorovas 
theory of local isotropy 

u^(d) = 2(e d)^/^ (3«1.4) 

By integration of the universal inertial subrange of the 
energy spectrum, it can be seen that only the effect of 
eddies whose scale is less than d is taken into account 
in the evaluation of u^(d) as given by (3,1,4). 

Hence the probability density of the relative 
velocity u between two points separated by distance d is 
given by, 

P(x) = — i exp [ 4l (5.1.5) 

^2% a 2 

where = u^(d). Hence the transitional breakage probability 
r (v) is obtained as 

r(v) = XPr[u2 = X erfc 

= Xerfc(Va(|)^/^ v“^AS/2 e ^^5) (3.1,6) 

where a = 2(2^^^-l) 

P 

Similarly 

r(d) = Xerfc (fa d"5/6/2 i 1/3) 


(3.1.7) 
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¥e define a maximtun stable drop diameter d^ such, that 


Ifa d;5/6 

21^73 


3.5 


( 3 . 1 . 8 ) 


as erfc(p) 0, if p ^ 3.5. 
•we have 


2(21/3-1) „l/3 62/3 


From squaring equation (5.1.8) 


a 

P 


(g)-l/3 4-5/3 

"71273 


12.25 


Thus, 




where 


C = — 

(2^/5~1)6 


For a baffled agitated vessel |^33] , 

e = 0^ 

Substituting, 


3t 2. 5 2/3 = 0, 1,2 i4/3 


1-5/3 _ 


= C-,C 


2 t 4/3 




1^3 


s 

( 1 ^) = 0 )- 0»6 _ n - 0.6 


L ^ cr 

- 0.6 


= ^4 


where = (^ 1 *^ 3 ) 


(3.1.9) 


(5.1.10) 


The functional form of this expression is identical 
to the empirical relations obtained for average drop diameter 
of equilibrium distribution in batch agitated liquid-liquj-d 
dispersions [ 3 I, 32 , 39 ]. 
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Redefining the transitional prohahility in terms of 
dimensionless drop volime and drop diameter we get, 

Y(-~) = Xerfc{3.5 (3.1.11a) 

s s 

y( d ^ xerfc{3.5 (-4-)“^^^} (3.1.11b) 

s s 

¥e recognise that this is a universal relation for 

all liq.uid dispersions whose viscosity and density are not 

very different from those of the continuous phase. The 

behaviour of y with is shown in Figure 3.1. 

^s 

From Figure 3.2, we see that the transitional breakage 

probability Y(-~) can be approximated by a series of 

^s 

'power law' approximations of varying exponents for different 

vol-ume regions. Figure 3.2 gives the region of dimensionless 

droplet volume and the corresponding exponent of 'power 

^s 

law' approximation of transitional breakage probability. 

From the experimental data of Madden and McCoy [54] on 
cumulative volume distribution of dispersed phase in a batch 
agitated vessel, Ramkrishna [4] has found that a 'power law' 
approximation for the transitional breakage probability 
holds good. The exponent obtained with respect to droplet 
volume V, by a preliminary analysis, was 2.0. This value 
agrees fairly well with the exponent 2.12 for the dimension- 
less droplet volume — range of 2 to 4.5 as can be seen 

^s 

from Figure 3.2. 





(l/A)/{v/vs) 




40 


Por the ’power law’ approximation of the transitional 
breakage probability, we have, 

r (v) = K v^ 

where K is the rate constant. 

Expressing the transitional breakage probability in 
terms of dimensionless droplet volume ( ^) > we have, 

Y( I-) = r(v) = b( 

. s s . 

where b is a universal constant for the region of dimension- 
less droplet volumes for which n is the value of the exponent. 
E = b v““ 

i.e. K a v~^ a d“^^ 
s s 


Prom the expression for maximum stable drop diameter d, 


we see that 


- 6 ^ ^ 
K a (s)^ ( P) ^ 


9n 


( 3 . 1 . 12 ) 


Prom the above expression for functional dependence of the 
rate constant K on P , e and cr, we see that K increases as 
e increases, and decreases as a increases. This behaviour 
is of course as it should be. 

Substituting for power input per -unit mass e in terms 
of agitator speed and impeller diameter in equation (3.1.12), 
wo get, 

18n 12n 9n - 9n 

KaII5 L5 p5 a 5 


(3.1.13) 
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Near d , the asymptotic ey-n.v, • 

Pansion of the complimentary 

error fmction yields an. a-D-D-nn^-: 

P oximate expression for y as 

X 1 


Y 


( A ') _ .. A.... . 1 / d X ^ /6 - ^ f’x 

^ HIT 3.5 ^ dj exp [-12.25 (-1-)"^'^^] 

^s 


For 


d 

d 

y(- 


(3.1.14) 


1 


s 


d 


^ ' F 3.5 ^ a^) exp[. 12 . 25 ( a )-5/3j ( 3 . 1.15 

- S 

d denoting an average diameter 

From Figure 3,3 » H can be seen +u -u ■ . 

^ hat this approximation is 

q.uite good upto g— = 1 , 5 . 

s shall later use this expressior 

of Y as given by equation (■31 nr-N 

• b; in the derivation of 
eq.uilibri-um drop volume distr-ih +• 

h ion in a batch agitated 

lean liquid-liquid dispersion 
Substituting for d in 

, , ® ®<luatlon ( 3 .I.U) from equation 

(3.1.10), the approximate ' 

"^ression for the transitional 
breakage probability near d n -u 

^ 'Can be written as, 

r(v) O 3 . If l 2/3 pi/a 5 /]_g - 

V exp[- n h 1 

II ^ 579 ^ 5^3 p J 

n ^ n -u • ^ ■ (3.1.16) 

and Cjj being constants. 

At this point, we compare 

present model for the 

transitional breakage probabii-i + 

- ■the model proposed by 

Coulaloglou and Tavlarides 

J* Their final expression for 

the transitional breakage 

■s ProbabUity is giy^^ 



/(d/ds)/^ 



Fig. 3-3 -Approximation tor transitional breakage probobil 
K near ds • 



43 


C 0 

T (T) = Cj y-2/9 N exp [ - ^ ' 

The fimctional form of this expression does not agree 

with the present model. However, comparing the above 

expression with the asymptotic expansion of transitional 

breakage probability near v^ as given by equation (3.1el6), 

we find that the functional forms are similar and that the 

arguements of the exponential function in both the cases are 

identical, Coulaloglou and Tavlarides [55] have not defined 

a maximum stable drop diameter. However, their functional 

dependence of maxim-um stable drop diameter on agitation speed 

H, viz, d a agrees with equation (5.1*10), They 

s 

assumed that the fraction of drops breaking is proportional 
to the fraction of turbulent eddies colliding with the 
droplet that have a turbulent kinetic energy greater than the 
droplet interfacial energy. They also assumed that the 
probability distribution of kinetic energy of eddies is 
exponential. Such an assumption is not supported by the 
experimental measurement of probability density of velocity 
at a point in agitated vessels [21, 23], They introduced 
the concept of 'breakage time' in their derivation of 
transitional breakage probability. The breakage time was 
estimated by assuming that the motion of the centres of 
mass of the daughter droplets to be formed by binary breakage 
is similar to the relative motion of two lumps of fluid in 
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a turbulent flow field. In their estimation of breakage 
time, they did not take into account the resistance to 
deformation and breakage of droplets because of the 
restoring force due to the interfacial tension. 

5.2 EQUILIBRIUM DROP VOLME DISTRIBUTION EOR BATCH AOITATED 
LEAN LIQUID-LIQUID DISPERSIONS; 

In an agitated batch liquid -liquid dispersion, 
continuous breakup and coalescence of drops occur. Eventu.ally, 
a dsmamic equilibrium between breakup and coalescence is 
established and a spectrum of drop sizes results. For a lean 
liquid-liquid dispersion, the rate of coalescence, though 
small, becomes important near equilibrium as the equilibriim 
is dynamic. Hence at equilibrium the rate of breakage, in 
order to compensate for the rate of coalescence, has to be 
small. So, we can expect the equilibrium distribution to be 
narrow around the maximum stable drop diameter d . Experi- 
mental observations of equilibrium distribution do show that 
they are narrow [30], 

A population balance equation for a batch liquid -liquid 
dispersion is given by, 

= I v(vO r(v*) p(t,v*) 0(v*,t) dv» - r(v) 0(v,t) 

V 

T 

^ 0(v’,t) iZS(v-v’ ,t) dv’ 

o 

00 

p 0(v,t) 0(v*,t) dv* , 


( 3 . 2 . 1 ) 
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whore 0(v,t) = .lumber density of drop population, 

v('v’) = Mean number of daughter droplets obtained 

by breakage of a parent droplet of volume v' 
p(v, v' )dv = Probability that a daughter droplet arising 

from a parent of volume v’ has its volume 
between v and v+dv 

and [i = Coalescence frequency. 

At equilibrium , 

GO 

0 = j v(v') r(v') p(v,v') 0(v') dv’ - r(v) 0{v) 

V 

V °° 

+5 / p 0(v’) 0(v~v') dv' - / p 0(v) 0(v') dv' 

° ‘ ° ( 3 . 2 . 2 ) 

¥e assume that the coalescence frequency p is a 
constant as the equilibrium distribution is narrow. Moreover, 
at equilibrium, the rate of breakage is small i.e. the energy 
of oscillation available from the turbulence of the continuous 
phase is only slightly hi^er than the minimum increase in 
the interfacial energy for breakage. Hence we can expect 
the breakage to be mostly binary. ¥e also assimae that the 
daughter droplet size distribution is uniform (i.e. perfectly 
random), 

i.e. v(v) = 2; p(v,v') = -^ 

Hence the popiilation balance equation is. 
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QD 

0 = 2 J r(v») -- r(v) 0(v) 

Y 

V 

■^2 0(v-‘V'O dv* - ji 0(v) 

o 


CO 

J 0(YOdv' 
o 

( 3 - 2 . 3 ) 


Bajpai, Eamlcrishiia and Prokop [47] have proposed an eq.uilibrinm 
drop size model based on viewing the coalescence-redispersion 
process as a single step process based on a constant coalescence 
frequency and a perfectly random redispersion. This model, 
which showed good agreement with experimental drop size 
distributions, is not applicable for very low dispersed phase 
fractions in which both breakage and coalescence rates are 
very small. This is the situation, which the present model 
is aimed at. 


As 0(v) = 




V 


, the population balance equation can 


be written in terms of volume density as, 

00 

f^^(v' ) r (v) f^(v) 


0 


2 J r (v» ) 


V 


,2 


dv’ - 


V 


+ 


li 

2 


/ 




f (v-v‘) f (v)°^f (v’) 

■ - F W ^ -V-/ 


(3.2.4) 

Let V be the volume at which the equilibrium density f-^C'vj 
exhibits a maximum. We define a non-dimensional drop 
volume X = — ^ . In terms of non-dimensional volume x 

V 

equation (5.2,4) is 
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For a narrow equilibri-um distribution, it can bo shown that 
(see Appendix l) the coalescence term 


ii 

2 



tlLs:3:i 

(5 


] dC 


! 


is of a much smaller order of magnitude than the other terms 
in the above eq,uation. Hence equation (3.2,7) becomes, 


0 = ^ ILK . .LU3.) -Y t( ? ) f(C ) - Y ( ? ) f'( 5) 

-p f(c ) K + ti ^ < 

® f(E ) 

where < = I !■ >. .. / g j _3 ^ constant. 

o 5^ 

Thus 



5b- K ~ Y ( ^ Y’ ( K ) 

C l\i- K +y{ ^ ) 


( 3 . 2 . 8 ) 


(3.2.9) 


In Appendix I, it is shown that b<2^Y(?) so that 
equation (3.2.9) becomes. 


|J(i) _ 5 \x:K-y(K) - g r'(g ) 

f(TT Cb fc 


(3.2.10) 


Since we expect the distribution to be narrow around 
d , we use the approximate expression for y near d given by 
equation (3.1,15) 

Y'( ^ ) = A exp [-a B (3.2.11) 

where A = (|“)^/^; B = 12.25; a =(|l)“5/3 ( 3 . 2 . 12 ) 

s s 

Substituting for y in equation (3.2.10), 
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^ exp [-a B A a B ^ 5 ~®/^exp[ -a B 


f-tH = -^ — ^ — — — ^ — 

(5»2.15) 

Since f ' (l) = 0 

5 IJ. < - A exp,[-a B] - A a B ^ ■exp[-a B] =0 (3 o 2.14) 

The elimination of A from equation, (3 •2.13) using (3,2.14) 
yields, 


= 5 [ C'^-exp [- aB( ?*'^/^-l) 


.-1 


*• 5/3 


Now :2::^ 18 and for C 1 we hare 

^ cc B E 4 . E / 

5 ^ ^ ^ ^-8/3- 

^ a B + 1 
3 


(|a B ? "®/^+ ?“^) 


_V 3 

(| a B + 1) 


-] 


(3.2.15) 


Hence 

= 5[ exp{ -a BU l) > ] (3-2.16) 

Integrating, 

In f( ^ ) = 5[ln5 - ^ exp {-a B( ^“^/^-l)} ] + C 

(3.2.17) 


where C is the constant of integration, E 2 panding In E and 
exp {-a B( ^~^/^-l)} in a Taylor series around E =1 and 
retaining terms up to second order, we obtain 

In f(5 ) = 5 [ - I (a B-1) (5 -l)2]+0 

Normalizing the non-dim ensionalized volume density and 
denoting it by f^( ^ ) , we have. 


(3.2.18) 
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^ 2 

^ [- " ] (3,2.19) 

^ a 2 

=25T^^) 

From equation (5.2.19), we find that the non-dimensional 
equilihrinm volume distribution is normal. Interestingly 
enough the empirical equilibrium drop volume density in 
terms of non-dim ensionali zed drop diameter obtained by 

Chen and Middleman [50] is also normal. Quantitative 
comparison of the derived distribution with the experimental 
data is possible only if the value of a is known. 

3.3 EVAIUATIOM OF ‘POWER LAW' COMSTAMTS OF TRAMSITIOMAL 

BREAKAGE PROBABILITY AND DAUGHTER DROPLET DISTRIBUTIOU 
FROM THE TRAHSIEHTS OF DROP VOLUME DISTRIBUTIOH IH BATCH 
AGITATED LEAH LIQUID-LIQUID DISPERSIOHj 

5,5,1 Equation for the Transients of Drop Volme Distribution 
in Batch Agitated Lean Liquid -Liquid Dispersions; 

We characterize the droplet population in a batch agitated 

lean liquid-liquid dispersion by the cumulative distribution 

function of drop volume F(v,t) defined as, 

F(v,t) = volume fraction of droplets whose volume is 
^ V at time t. 

For a batch agitated lean liquid-liquid dispersion, it can 
be seen [4] that, 

= f r(v') G(v,v') dP(v‘,t) 

V 


(5,3.1) 
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In eq.nation (3.5.1), the right hand side is a stieljes 

integral which may also be written as, 

CD 

r (v») &(v,v«) (iv' (3o3o2) 

where is the density fnnction. 

Eq,uation (3.3.1) may be solved for E(v,t) if T (v) 
and G-(v,v’ ) are known. However, our present interest is one 
of determining r (v) and G(v,v') when F(v,t) is experimentally 
measured. 

If the 'power law' approximation for the transitional 
breakage probability is valid, then 

r (v) = K v^ (3.3.3) 

For the range of drop sizes encountered in agitated 
liq.uid-liq.uid dispersions, breakage is not Chaotic [8] and 
is caused by oscillations of droplets due to random velocity 
and pressure fluctuations [8]. It is believed that the 
mechanism of breakage is the same for all droplets in which 
case one would expect the daughter droplet distribution 
Gr(v,v') to be a function only of the ratio of daughter and 
parent droplet volumes ;i.e, once a droplet is broken 
the daughter droplet distribution is independent of parent 
droplet volume. We call this 'similar breakage'. If the 
breakage is 'similar*, then 

G(v,v') = g( 


(3c3.4) 
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Then equation (3»3.1) hecomes, 

QO 

• = J K g( |,) dF(v»t) (3.3.5) 

9 V 

3.3.2 Estimation of ’Power Law' Constant n; 

Fillipov [53] has shown that the solution of equation 
(3.3.5) under the assumptions (3.3.3) and (3.3.4) e-venfcually 
approaches an asymptotic distribution in the similarity 
variable s = (1+Et)v^ That the grouping (1+Kt)v“’ is 
indeed a similarity variable is easily verified by letting 
f(z) = F(v,t) and substituting into equation (3.3.5) which 
yields, ^ 

Zf'(2)= / Cf’ ( C ) gn (f-) dc (3.3.6) 

where f is the derivative of f and 

C = (1+ Et) v'^ ; = g[(-|~)^/'^]. 

Equation (5.3.6) can be rewritten as, 

CO 

0(z) = / (!(( 5 ) s, (^) dc (3.5.7) 

where 0(z) = z f ’ ( z ) 

For times sufficiently large, the similarity variable 
z is nearly equal to Etv^ which suggests an interesting 
possibility of evaluating the 'power law' constants E and n 
from the measurements of cumulative volume distribution 
F(v,t) in batch systems [4]. If for large times, Eiv*^ 



53 


is indeed the similarity transformation, the cumulative 
drop volume distribution P(v,t) for different times when 
plotted in terms of the similarity variable should collapse 
into a single distribution. Therefore, for a fixed P(v,t), 
the corresponding values of v and t shou.ld be such that 

Etv^ = constantj i. e. v^ t = constant. 

Thus a log-log plot of v versus t for a fixed ^(Vyt) shouid 
yield a straight line with a slope of - Further, similar 
plots for different fixed values of F(v,t) should yield a 
set of parallel lines if the assumptions (3.3^3) and (3.3,4) 
are tenable, A preliminary analysis of the transients of 
drop volume distribution obtained by Madden and McCoy [54] 
by Ramkrishna [4] does yield such parallel straight lines 
thus vindicating the assumptions, 

3*3«3 Evaluation of K and g(x); 

For sufficiontly large times, one may use the similarity 
variable y=v^t which gives, 

CD 

y f*(y) = Kj^fU? ) gT(f^)(iC (3.3o6a) 

y -1- 4 

Multiplying (3 *3. 6a) by y^”^ and integrating from 0 to “ 
one gets, 

oo CD 

/ y^ f’(y) dy = Ef y^"^ / Cf’( ?) gT(-^)d^ dy 

° ° CD y e ^ 1 ^ 

== K J^f'C c) J y gT(-^)<iy (3o5o8) 

0 0 ? 
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i.e. - K 

CO _ 

4^^ = J y f ' (y) <iy 


where 


1 T. 

and = J X g^(x) 


dx 


(3e3.9) 

(3o3.10) 

(3o3«ll) 


variable y 


If we define a function in terms of a new similarity 
= V t^/^ such that 


(y^^^) = f (y) 


with moments 


00 


if (y^/“) 

O 

then it follows that 
^t/xi ~ ^r 


(3c3.12) 

(3.3ol3) 

(3.3.14) 


If we define ^ 

= I ^ (3.3.15) 

0 

where { ) a^e the moments of the daughter droplet distri- 

bution g(x) y 

then, it is easily seen that 

^ Y(r+l)n-l (3.3.16) 

low, equation (3.3.9) may be transformed to yield, 

^ ^r+n P(| -l) (3.3.17) 

which in combination with (3.3.16) gives, 

= K n X^_^^ Yj._i 


(3.3.18) 



55 


Measurements of tiie transients of the emulative 
volume distribution F(v,t) in a batch system shoTuLd in 
principle give the moments { } after transformation of the 

distribution using the similarity transformation y^'^^ = v 
Thus, from a plot of versus y^-^^ the momonts[ p., } 

can bo obtained. 

Equation (5.5.18) then provides the values of {K )■ • 
From the values of {K i't is possible to estimate the 
function K g(x) by means of a suitable orthonormal expansion 
[43]. Since g(l) = 1, K can be evaluated from the value of 
K g(x) at x=l. However, from equation (3.3.18) it can be 
seen that the evaluation of moments { y^} involves knowledge 
of moments of d Y /dy^/^ of higher order. Due to inaccuracies 
involved in the evaluation of higher moments of experimentally 
measured distribution F(v,t), it may not be possible to 
evaluate the higher moments of g(x) accurately. However, 
since g(x) ranges between 0 and 1, the moments of 
progressively decrease. Hence, it may be sufficient to 
evaluate only the leading moments of g(x) in order to evaluate 
g(x) satisfactorily throu^ polynomial expansion. Further- 
more, zeroth moment of g(x) represents the average volume 
of the daughter droplet and the first moment the standard 
deviation around the mean. This in itself constitutes useful 
information about the daughter droplet distribution. 



The f motion K g(s:) can be expressed as a sm of 
orthonormal polynomials, the coefficients of which are 
related to the moments {h y^} . The orthonormal polynomials 
can be generated with respect to a weighting fmction w(x) 
which is a good approximation of g(x) i.e. orthonormal 
polynomials 0^(x) can be generated by orthonormalising 
1, X, x^ , , . , with respect to the weighting fmction w(x) 
such that, 

1 

/ w(x) 0^(x) 0j^(x) dx = 0 n 0 m 

° =1 n = m 


Expanding K g(x) in terms of {0^(x)} , we have. 


GO 


K g(x) = I w(x) 0^(x) 


then, 


n=0 

CO 


n' 


(3.5.19) 


= 4 K Y^, K Y^ K (5.3.20) 

Thus, the polynomial expansion introduces a series of 
perturbations aromd the weighting function w(x) the 
coefficients of the perturbations being deteimined by the 
moments { K y^ } [45]. If the weighting fmction is a good 
approximation of g(x) , it may be sufficient to consider only 
a few terms in the expansion i.e, the expansion will converge 
after the inclusion of a few terms. 



CHAPTER 4 


A STATISTICAL THEEMODYNMIC MODEL EOR EQUILIBRIUM DROP SIZE 
DISTRIBUTION IN AGITATED LEAN LIQUID-LIQUID DISPERSIONS 

Dispersed phase systems are of common occurrence in 
many industrially important processes. The study of agitated 
liquid-liquid dispersions has heen the subject of nimerous 
investigations. Due to droplet phenomena of breakage and 
coalescence, dispersed systems reach a d 3 ?namic equilibrium 
wherein the drop size distribution is invariant with time* 
Knowledge of the equilibrium drop size distribution is 
important for characterisation of such systems, 

A statistical mechanical formulation for such systems 
was proposed by Hulbert and Katz [ 43 ] where the system was 
described by a population balance equation. Population balance 
equations were used by Valentes et al, [3, 46] to describe 
the behaviour of agitated liquid-liquid dispersions. Breakage 
and coalescence, being random due to random velocity and 
pressirre fluctuations in the tixrbulent flow field, were 
represented by related probabilities, namely, the transitional 
breakage probability, coalescence frequency and size distri- 
butions for newly formed droplets by breakage. The prediction 
of equilibrium drop size distribution using such a kinetic 
approach necessitates knowledge of breakage and coalescence. 
Even for lean dispersions, where the interaction between 
droplets is negligible, it is necessary to take coalescence 
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into account for predicring the eq,uilibriuiii drop size 
distribution from a kinetic approach. This makes the 
population balance equation non-linear due to which it becomes 
necessary to make simplifying assunptions for the equation 
to be amenable to solution. Curl [ 44 ] and Erickson [57] 
assuned that the droplets redisperse immediately after 
coalescence without change in size. Bajpai, Ramkrishna and 
Prokop [47] proposed a coalescence-redispersion model for 
drop size distribution in agitated, vessels by assming that 
droplets coalesce in pairs and redisperse instantly into new 
pairs arbitarily. In Chapter 3, the equilibriua drop size 
distribution in an agitated lean liquid -liquid dispersion 
was predicted using population balance equation under the 
assumptions of binary breakage, perfectly random sizes for 
daughter droplets from breakage and a constant coalescence 
frequency. Although, justifications for these assumptions 
can be given, the prediction of equilibrium drop size distri- 
bution hinges critically on them. 

In this chapter we investigate the problem of predicting 
equilibrium drop size distributions from a statistical themo- 
dynamic viewpoint which considers only the total energy of the 
system i.e, the kinetic, potential and interfacial energies 
of the droplets. Since the thermodynamic approach does not 
explicitly take breakage and coalescence into account, it is 
not necessary to make simplifying assixaptions in regard to them. 
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Moreover, in case of dilute dispersions the thermodynamic 
approach neglects the effect of interaction between droplets 
as it is negligible. This is justified as the phase average 
of th.e energy of the system is a very weak function of 
interactions [58], though droplet interaction is, in principle, 
very important without -viiich the system would never attain 
a dynamic equilibrium, 

4.1 THE SYSTEM: 

The system consists of the dispersed phase in the form 
of droplets of various sizes exposed to a turbulent flow 
field in a well-stirred vessel. The flow field is stationary 
i,e, the droplets are exposed to a turbulent flow field 
whose characteristics are time invariant. The flow field 
can also be considered as homogeneous and locally isotropic 
[20, 21, 23, 24, 28], Eventhou^ the droplets continually 
interact with their environment i,e, the continuous phase, 
we would consider the system to be consisting of only the 
dispersed phase in the form of droplets. This is justified 
as we will see later that the interaction of droplets with the 
continuous phase is accounted for by considering (i) the 
convective motion of the droplets (ii) the oscillations of 
the droplets and (iii) the potential energy of the droplet 

^Eventhough inhomogeneity is observed in a well-stirred vessel 
due to most of the energy dissipation taking place in the 
impeller stream [18], the flow field is fairly homogeneous 
outside the impeller stream. 
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at its eq.uili’bri'uiji spherical shape. The power input into 
the system per unit mass e, total volume of the system 
V(volume of the dispersion) and dispersed phase fraction 0 
are fixed. Since the total number of droplets is determined 
by breakage and coalescence, one would expect the total 
number of droplets for fixed e, V and 0 to be random as 
breakage and coalescence are random. Since the microscopic 
system consists of a very large number of droplets, one 
would also expect the fluctuation in the total number of 
droplets around the average (ensemble average) number U to 
be extremely small. In other words, one can, for all practical 
purposes, consider the total number of droplets in the system 
at equilibrium to be constant. As the interaction between 
droplets is negligible for a lean liquid-liquid dispersion, 
the system can be assmed to consist of N independent droplets. 
This is justified as the phase average of the energy of the 
system is a very weak function of interaction [58]. In 
order to derive the equilibrim drop size distribution, we 

( 

consider a canonical ensemble of subsystems, each subsystem 
being a droplet i.e, each droplet is exposed to a constant 
'temperature*. 'Temperature' here refers to the characteristic 
of the turbulent flow field. Therefore, we need only to ' 

evaluate the partition function of a droplet in order to 
evaluate the equilibrium drop size distribution. 
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4,2 PARTITION FUNCTION OF A DROPLET; 

A droplet exposed to a turlDulent flow field will be 
acted upon by turbulent eddies of different scales or 
frequencies, Eddies of scale greater than the drop diameter 
will impart a pulsating convective motion to the droplet. 

There will be relative velocity and pressure fluctuations on 
the surface of the droplet due to eddies of scale less 
than the diameter of the droplet which will result in oscillatio 
of the droplet around its spherical equilibrium shape. Because 
of the interface, the droplet will also have potential energy 
which will be a function of the interfacial tension and the 
interfacial area. 

Therefore, the energy of a droplet can be split up 

into : 

(i) Kinetic energy of pulsating convective motion 
(ii) Kinetic and potential energies of oscillation 
and (iii) Potential energy at its spherical equilibrium 
shape. 

4.2,1 Kinetic Energy of Convective Motion ; 

The convective piilsating motion of a droplet of 
size X is caused by eddies of scale larger than X. Since 
the drop sizes lie mainly in the inertial subrange of the 
universal energy spectrum [ 28] , the total t\jrbulent energy 
available for pulsating motion exhibits very small variations 
in the prevailing size range. Therefore, the convective 
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pulsating motion can be assmed to be independent of drop 
size [48], Therefore, 

Kinetic Energy = (p^ + Py + P^) (4*2.1) 

where p^, Py, p^ are the momenta corresponding to the 
cartesian coordinates x, y and z and M is the mass of the 
droplet, 

4.2,2 Energy of Oscillation ; 

A spherical droplet is acted upon by turbulent 
velocity and pressure fluctuations on its surface. The turbu- 
lent velocity and pressure fluctuations can be considered 
to be due to eddies ‘striking* the surface of the droplet. 

The eddies (or velocity and pressure fluctuations) create 
oscillations of the droplet and hence their defoimation. Hear 
equilibrium, the rate of breakage is small i.e, the defoimation 
of the droplets is snail. This is because the energy of eddies 
of scale smaller than the droplet diameter is of the same 
order of magnitude as the minimum energy required for breakage. 
Moreover, small size droplets more or less behave as rigid 
spheres i.e. the tangential momentum transfer across the 
interface resulting in internal circulation is small. There- 
fore, the droplets can be assumed to undergo pure deformation 

without rotation! Hence the induced flow field inside the 

o' ' ' ' ■ 

This assumption is only approximate as there will be some 
internal circulation. 
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droplet due to its oscillation can be considered to be 
irrotational. For small pure defoimation, a sphere will be 
deformed into an ellipsoid. It is evident that for oscillations 
of a droplet, the velocity and pressure fluctuation can occur 
only at a point on the surface of a droplet at a time, as 
otherwise, the effect of velocity and pressure fluctuations 
would be to create bulgy deformation of the droplet [8]. Hence,- 
the oscillations of a droplet are to be viewed as forced 
oscillations of a sphere due to a random periodic disturbing 
pressure acting at two diametrically opposite points on the 
surface. Since the equilibriim drop size distribution is 
narrow and the drop sizes lie in the inertial subrange of 
the universal energy spectrum, we can assume the turbulent 
energy available for oscillation to be independent of drop 
sizes. 

Consider a droplet of size X, At any time, the droplet 
will be acted upon by a random periodic pressure of random 



frequency at any two diametrically 
opposite points A and B. Since 
the droplet is incompressible, the 
effect of this pressure will be to 
induce a pressure distribution 
p(9, 0) along the surface of the 


droplet. Because of local isotropy, the droplet will oscillate 


between oblate and prolate spheroidal shapes. Therefore, it 
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should he possible to find out the pressure distribution 
P(&, 0). 

The periodic pressure fluctuation p(&, 0) can be 
written as, 

P (©, 0) = Pa 0 ) e~^ 
w being the frequency of the fluctuation. 

Here both p„(©, 0) and w are random variables. Since the 
deformation of the droplet is small, the function 0) 

can be represented in terms of surface harmonic functions as, 

Pa(®» 0) = I I cos n0 + sin n0) (cos ©) 

m n 

^ ( 4 . 2 . 2 ) 

where P” (cos ©) = sin © , is what is 

"" d (cos 9)"^ 

called an associated Legendre function, P^(cos 9) being the 
Legendre polynomial of order m. As is well known, m takes 
all positive integer values from zero, while n takes the 
values 0, +1, +2, , . , + m. 

The constants and are given by, 

2-11; m 

“mo “ 1 1 e ae d0 ■ (4.2.3) 

% 

f P.(e,«() { cos 9 ) cos n0sin9 dO d 

Q in 

(4.2.4) i 

% 

i P„(©»0) I^(cos 9 ) sin n0 sin 9 d9 dj 
0 ^ 

(4.2.5) 




- 2 m+l (m~n) I 
mn 2 % (m+n) I 


2 % 


0 


and 

^mn 


2m +1 
2 % 


(m~h) I 
(m+n) i 


2n 

! 

0 
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Once Pg^(©, 0) is known, it is possiTole to evaluate the 

constants a _ and B . 

mn '^mn 

The equilibrium spherical surface of the droplet is 
given by r=sR (droplet radius). We can describe the deformed 
surface of the droplet by the equation r=R + 5(^,0), where 

C(&>0) is the deformation of the droplet from its spherical 
shape. As stated earlier, the induced flow field due to the 
oscillation of the droplet is assumed to be irrotational and 
hence can be described by a velocity potential Y'(r, ©, 0). 

The velocity potential Y*(r, ©, 0) satisfies the 
Laplace equation 

y ~ 0 (4.2.6) 

with the boiindary condition [59], 


9 '?' ,2 

‘’Tt + 


2L 

2 


R 


1 

T 






+ Pq + p(©,0) = 0 


(4.2.7) 


at r = R 

where p^ is the average pressure in the continuous phase. 

Defining a new velocity potential + ^o^^* 

we get, 

- 0 (4.2.6a) 


with the boundary condition. 



[25 + { 


_1 

sin© 


is (sine |S^) + 




sin © 3 0' 


> 1 


+ P(», 0) = 0 


at r=R 


(4.2.7a) 
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Differentiating (4.2.7a) with, respect to time and putting 

3 W 3 V 

— ^ , we have the modified "boundary condition. 


3^ y g r p _3I. . X- f -1 3 

^gt^” ar ar^ sin© © 

+ P(9, 0) = 0 at r=R 


(sin© -||) 


jil 

sin^© a 0^ 

(4,2.7b) 


}] 


How, 

p(©, 0) = w Pg^(©, 0) 

Since ¥ satisfies Laplace eq.uation, 'F can be represented as 
a linear combination of volume spherical harmonic functions 
sin ^ (cos ©), Therefore we seek a solution of 
the foim, ' 

'^ = I I 3:*“ (a^^^ cos-n0 + b^^^ sin n0) (cos©) 

(4.2,8) 

Substituting for ¥ in equation (4.2.7b), we get, ' 

p w^ V2) ^ j j ^ ^ 0 ^ ^£(cos ©’ 

m n 

^ I I (^mn ■*■ \n ^ 

m n i 

»(m-l) (m+2)] ”^2^ 

+w J I(%jj oos nJ* + sin nj!) :^(oos0)]=O; 

m n 


(4.2.9) 
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after using the relation [59] satisfied by the surface 
harmonic functions 


1 

sin © 
where 


3 

d© 


3 Y 

(sin© + 


3 © 
cos n0 
mn sin n0 


Sin^© 


32 V 


3 0' 


+ n(m+l) Y^^ = 0 
' mn 


(cos ©). 


( 4 . 2 . 10 ) 


m 


Therefore, eq.uation (4.2,9) reduces to 

I I { - P w2 ■-2j-E“m(m~l)(m+2)aj^^ 


+ w } + sin n0 { - p w^ + -% m(m-l) (m+2)b^^ 

] = 0 (4,2.11) 

Therefore , 

o = E «. 

mn m ^ O rr 

[ p w"^ ^ m(m-l) (m+2)] 

R^ 


w 


PR“^[w^ - w2 ] 
L mn ^ 


cc^ 


mn 


( 4 . 2 . 12 ) 


where w = 
nm 


PR' 


m(m-l) (m+2), w^ being the mth natural 


frequency of oscillation of the droplet [59]. 
Similarly, 


w 


mn opm r.2 „2 t 

PR [w - w^] 


mn 


Obviously, the relations (4.2.12) and 


(4.2.15) 
,2.13) hold good only 
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i'J 1 : ^ ¥ m. Measurement of equilibrium drop size 

d is trio lit ion for uater - CCl^ + i-octane (50-50 per cent) 

33^8 tem has shown that t3rpical droplet volume s:r lO"" cc. 

As cr c; 40 dynes/cm, the smallest natural frequency of 
oscillation w^2 ^ ^ 10“^ . Measurement of the energy spectrum 

in stirred vessels [25] shows that the inertial sub-range of 
the universal energy spectrum (-5/3 slope) corresponds to 
eddies whose frequency is of the order of 500 Hertz, There- 
fore, the natural frequency of oscillation of droplets is 
much higher than the forcing frequency of the random pressure 
fluctuation. 

Therefore, the velocity potential f can be expressed 
as 

'i' = I ^^^mn ^mn ^ (cos©) ■** 

m n 

From equation (4.2,3), we get, 

2b % 

“oo = iV i ® (4.2.14) 

The integral' is the total radial force acting on the droplet 
due to the induced pressure distribution p„(©, 0). The total 
radial force would result in radial oscillations of the 
liquid droplet. Since the droplet is incompressible, radial 
oscillations are not possible and hence the total radial force 
has to be zero. i,e. = 0. 
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Prom equations ( 4 ., 2,4) and (4,2.5), we get. 


“11 =5^ 


2 % 

71 

1 

21 

J 

JPa (^»0) 5j(cos ©) COS0 sin© d© d0 


0 

0 

Pll = 2^ 


2 % 

% 

1 

2 r 

J 

^ Pg^(©',0) (cos©) sin0 sin© d© d0 


271 

% 


“10 = 2^ 

J 

r\ 

f 

r\ 

Pa(©,0) (cos©) sin© d© d0 


and 

^10 “ 

Since cos©) = i sin© and (cos©) » cos©, 

2% Tl 

“11 ~ (4.2,15) 

0 0 
2tc n 

“10 ~ 2^ ^ ^a (4.2,16) 


and 


2 % % 

^11 “ " ' ^ “ " f f Pg^(©»0) sin© sin0 sin© d© d0 (4.2,17) 
o o 


It can be seen that the three integrals are the 
components of the total force acting on the droplet along 
X, z and y coordinates respectively. These components would 
result in the translation of the droplet. Since the induced 
pressure distribution p(©,0) results only in the oscillations 
of the droplet, the three components of the force are zero, 

i.e, = 0 . 



70 


Therefore, the summation in the expression for ^ 
starts only from m=2 i.e, m=2 corresponds to the smallest 
frequency of oscillation. 

The kinetic energy of oscillation T is given hy, 

^ = ///(V'i'-V¥)dv 

V 

By applying Green’s formula, we get, 



where n denotes the outward normal to the surface. 

Therefore , 

2 % % I 

T = |p / / j sin© d9 d0 (4.2,18) 

o o r=S 


substituting for 'P from equation (4.2,8) and using the 
following relations, 

2n n 

/ / ^pq ° 

o o 

2n n 

/ / (Q»0) (^,0) sin© d© d0 = 0 m,n 0 p,q. 

0 0 


and 


2m 71 _ 

^ f C&»0)]^ sin© d© d0 

0 0 


/ 4 It 
2n+l 

2m 

'-’2n+l 


m=0 , i=e 
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(cos9) cos m0 

^mn ^ (cos©) sin m0 

and i=e or o , we get. 



I [m E 

m=2 


2m+l 


4 % 
2ni+l 


a2 + 
mo 


I ( 

n=+l.. , 


(a^ + )) } 1 

mn mn^' •* 


^-2i(wt + 7 i/2) 


2n (m+n) I 
2m+l (ia-n) I 


( 4 . 2 . 19 ) 


We can define the momentum of oscillation 


•Ti 


by the relation, 


1 

2 M 



= T. 


(4.2.20) 


Therefore, 

= [M P i [m E.2m+1 
m=2 


2m+l ^mo 


I 

n=+l 


/ 2 It ( m+n) 
'2m+l (m-n 


(a^ + b^ ))n^^^ 
'■ mn mn^ ' ^ 


( 4 . 2 . 21 ) 


The velocity of oscillation is given by, 

7=^21 
os M 


( 4 . 2 . 22 ) 


Since a^^ , b ^^ and w are random variables, and 7^^ are 
also random variables. 

Substituting for Y from equation (4.2.8) -in equation, 


(4.2.7a), we get. 
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m L ^mn (cos9) e”^ ^ 


m 


- -^ [2 5 + ^ (sin© |1^) + 


3?L 


] 


“ ^ sin^9 3 0^ 

\ ^^“mn + ^mn I^CcosG) e~^ ^ = 0 

(4.2. 23) 

¥e seek a solution for 5 of the form, 

^ = I + ^mn sin n0) ig(cose) e“^^ 

m n 

(4.2.24) 

Substituting for 5 in equation (4.2,23) and using the relation 
(4*2,10), we get, 

I + -75 (»«)(«-!) V %n > “f 


m n 


R 


+ { pw E (m+2)(m-l) sin n0] (cosO) 


= 0 

-a m 

i.e. c = -^mn 
mn 


—2 ' + -rS — 3 

^ ^nm ~^L>' 


(4.2.25) 


and 


*“Pr«v, 21 “ 2 

d_.. = . .. [ 1 + w 


^ E w2 

run 




(4.2.26) 


run' 


lake a surface element between 9, 9H-d9 and 0 and 
0 + d0. Ihe area of the surface element df is given by, 

df = f(E+5 )2 + I f (^)2 + _1 ( |_|)23 de d0. 

sin 9 ^ 
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The area of the spherical surface element df^ is given hy, 
hfg = R sin9 d& d0. Therefore, the increase in the area 
of the surface element due to deformation is given by, 

idf = df - df^ = [ZRK + 5^ + I { (-||)2 H- (-|f)2 ] 

BXH 

sin© d© d0. 

3 ? 3 5 

Since -g-g and are small compared to § , neglecting the 

terms within the braces on the right hand side, we get, 

Adf = (2R? + sin© d© d0. 


Ror small deformation, we can take the increase in 
the potential energy or equivalently, the potential energy 
of oscillation, as the increase in the interfacial energy 
due to deformation* Therefore, the potential energy of 
oscillation is given by, 

= o f / (2R 5 + 5^) sin© d© d0. 

0 0 

2n % 

J / 5 sin© d© d0 

0 o 


2 % % 

= f I 1 I ( cos n0 + sin n0) 
o o m n 

£j(cos©) sin© d© d0 

-ivrt -lirt ^ n 


= 4ir c 
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Therefore, 

2% n 

P^g = cr J I sin9 d© d0 
o o 

Substituting for ^ from (4.2,24), 


(4.2.27) 


" H ^ “mo + 


OS 


m 


4% r.2 ^ Y r 2 It (m+n)l (^2 ^2 

^ ^ •C2m+r) (m-n) I ^ mn'^mn' ^ ^ 


n=+l . . . 
-2iirt 


(4.2.28) 


Defining the generalised deformation coordinate r\ by the 
relation, 

P = i a (4.2.29) 

os 2 ' 


„ r r 8 TC f>2 4. V / , 4 7 C (mfU ) I (o^ 4. }} T 

^ i ( 2 r +1) ^mo i ( ^+1) (m-n)T ^ mn ^mn-' 


m 


n=+l , . , 


-iwt 


(4.2.30) 


Therefore the kinetic and potential energies of oscillation 

are given by, 

[P — ..JL.. p2 . p _ 4 O’ 

^ - 2 M ^x] * -^os 2 ^ ^ 

where p is the generalised defoimation coordinate and is 
the momentum conjugate to p and both p and are random. 


4.2,3 Potential Energy of a Spherical Droplet ; 

The potential energy of a droplet at its eciuilibrium 
spherical shape is only due to the interface. It is assumed 
that the density and viscosity of the dispersed phase are the 



75 


same as tiiat of the continuous phase. At equilibri'um , the 
temperature and chemical potential in the droplet will he 
equal to the temperature and chemical potential in the 
continuous phase. Since the potential energy of a spherical 
droplet is only due to the interface, it should be equal 
to the work done in transferring a "volume v(= ^ 7 ?) of the 
dispersed phase reversi bly from a bath- of dispersed phase 
liquid at the same temperature, pressure and chemical potential 
as that of the continuous phase into the continuous phase 
in the form of a spherical droplet of diameter X replacing 
an equal volume of continuous phase. 

The internal energy and work function of the 
dispersed phase liquid of volume v in the dispersed phase 
liq-uid bath are given by, 


= T.S - P^v + p n (4.2..^S1) 

A^ = ~P2_v + p n (4.2.32) 

where = continuous phase pressure 
n = number of molecules, 

Now, the internal energy E2 and work function A2 of the 
droplet of vol-ume v are given by, 


E^ = T,S - P2V + o A + p n 
A2 = ■“^2'^ + cr A + p n 
where P2 = Pressure inside the droplet 


A = Interfacial area. 


(4.2.33) 

(4.2.34) 
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Therefore, A 2 -A^ = AA = - (^’ 2 -^ 1 ) v + a A (4.2.35) 

“ aA = Maximm work obtainable from the system in an 

isothermal reversible change from state 1 to 2. 
Therefore, the potential energy of a spherical droplet of 1 
is given by, 

(pTe)^^ = AA = -(P2-P;L^ V + cr a 

= -AP|x^ + cr7tX^ (4.2.36) 

Since a droplet at its eq.\ailibri-um spherical shape does not 
have any kinetic energy, the total energy of the droplet is 
its potential energy. Therefore, the Hamiltonian for a 
spherical droplet of size X can be written as, 

H = -( AP) I X^ + c -n; X^ (4.2.37) 

where X is the generalized coordinate, the drop diameter. If 
Pj is the moment-urn conjugate to the generalized coordinate X, 
we know that, 


iii - ■?, 

3X ^X » 

Substituting for H, 


3 H 
3PX 


3 H _ 

3 X 


A P) I X^ + 2 It c X. 


At equilibrium, the pressure difference between the droplet 
and the continuous phase is the capillary pressure p^. 
Therefore, 


AP 



(4.2.38) 
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Substituting, 


9 H 

TT = 


~4 cf % ^2 


T- + 201 crZ = 0. 


X 2 

Therefore =0 i.e. = constant. 


And 


9 H 


= 0 


-.e. 1 = 0 


or ^ • 


\J \JJL. 

X 

i.e. the generalized momenta p^ conjugate to the generalized 
coordinate X is identically zero for all values of X, This 
is understandable as a droplet at its spherical equilibrium 
shape does not have any kinetic energy. Substituting for 
P from (4,2.38) into (4.2,37), we get, 

(4.2.39) 


1 _ _ ^2 


(i’.E)e^ = f cr It X 


4.2,4 Perivation of Partition Punction of a Droplet ; 

The coordinates of the phase ^ace are x, y, z, t), X, 
Px» Py» ^T) ^X* ^ droplet has 6 degrees of freedom. 

The total energy of a droplet of diameter X can be written as 


E = 


2 M 


/_2 , 2 , 1 „2 , o 2 , G -^2. 

(P;^ + Py + P^) + 2lf + 2 ^ X 


(4.2.40) 

Therefore, the Hamiltonian H of a droplet can be written as, 

^ - 1 <'-o2 . ^2\ 1 ^2 ^ O „2 , o ^ y2 

H(X, p ) - 2 M ^y '^2 M^ti‘^ 2^ + 37tX 

(4.2.41) 

Hence, the partition function of a droplet can be written as, 
QC7, ?)=//,.. / exp [-p (p^ + p| + Pf) + ^P^ 


12 M '-^x " -^y ‘ • 2M -^n 

+ I T]^ + I It X^ } ] dx dy dz dp^ dp^ dp^ dp dp^ dX 
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where ~co < p^, < co , -co h < go , -co < p^ < co , 

0< X < CO . 

Here the iiategration is not performed with respect to pj as 


p^- = 0 ¥ X. 

Integrating with respect to x, y and z, 

P) =7^ exp[~p { 2 "^ (Pjj + Py + Pg) 2 M ^ 

+ f + I n: } ] dp^ dPy dp^ dt) dp^ dX. 

How, M = oc X^ where a = ^. 


and 


-OD 


-00 


CD - *^ / 

/ exp [- — L p2^ dp^ = (^2L|JL)1/2 


2a X 


CO 


/ exp [- I cr -n^] dTi = ( 


OD 


therefore, 

Q(V,p) = V J exp [-p | a X^^dX 


= V 


15 37^ 1 


TC C) 


How, 


■ Xa. ..Q = 6 ^ g 


(4,2.42) 


where E = Average energy of a droplet. 

Therefore, p = J 
E 

If we define a ‘Dispersion Temperature’ T^ such that, 

1 E 


(4,2.43) 


1- « c ^ 


T = I 

■^d 


(4.2.44) 
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where k is the Bolt 2 mami’s constant. 

It is interesting to note that the average energy 
of a droplet E = 6k If the 6 degrees of freedom of 

a droplet could be further partitioned i.e. independent, 
then one would have expected the average energy of a droplet 
E to be 5 according to the law of eq.uipartition of energy. 

Erom equation (4,2,44) it is clear that the total energy of 
the droplet is not equally partitioned into different degrees 
of freedom. This is under standard able as the kinetic energies 
of convective motion and oscillation of a droplet are not 
only functions of the corresponding momenta but also functions 
of droplet diameter, a generalized coordinate i.e. the 
different forms of energy of a droplet are not independent. 


4,3 DEEITATION OE EQUILIBRIUM DROP VOLUME DISTRIBUTION; 

The probability density of drop size f(X) is given 

]■/.../ exp[-p H] dx dy dz dp^ djy dp^ du dp^ 

f (X) = ^ 

Q (V, p) 


Evaluating the numerator and substituting for Q(V, p), we get, 


f(X) = 


X® exp [-P 






(_6_)7/2 J/2 2^ x2;| 

E E 

(4,3.1) 
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Let Z* = 

'■2 Ti: a' 


(4.3.2) 


ITon-dimensionalizing f(Z) with respect to X and taking into 


account the relation f(-^) = f(X), we get, 


Z 


f(-^) = ^ 

Z 


_1 / Z x6 ( X ^2^ 

15 L/2 ^ exp[— ( ' ' jut') J 


(4.3.3) 


% X ■ - X 

The equilihrium drop volume density in terms of drop diameter 
f^ (Z) is given by. 


V 


(X) = .^( a . 

^ fxh(X) 


ITon-dimensionalizing f^(X) with respect to X 


r - X X . 


X A 2. 


(4.3.4) 


Now, 


O' 


CO 

/ X^f (X)dX 


J X‘^f(X)dX 


0^ 


Mi., 


X' 


“00 


J f(2.) a( 5.) 

Ox" X X 

T* 4^ 4- 1 . , 

= ^ 43 ^2 2 

46 4 

15 2 1 


^* 16.24 4_ 

^ 35.3 ^ 4/2 


i.e. d^2 = A X 


where A = -^2 


(4.3.5) 
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Therefore , 


d 




X 


'■32 


X 


^10 


X 


X* 


i.e. 


= ^ ( ^)9 exp [-a2(J^)2] 


( 4 . 3 . 6 ) 


'32 


'32 


'32 


Pigure 4.1 shows 

density — ) 

<132 

(=^). 

<3.32 


a plot of non-dimensionalized drop volume 
versus non-dimensionalized drop diameter 






CHAPTER 5 


EXPERIMENT Ali TECHNIQUES, APPARATUS AND PROCEDURE 

5.1 EXPERIMENTAL TECHNIQUE: 

The droplets in an agitated lean liciuid -liquid 
dispersion were encapsulated with a thin nylon film at the time 
of measurement through an interfacial poly condensation reaction. 
The encapsiuLation stops further breakage and coalescence of 
droplets thus freezing the population at the time of measure- 
ment, The encapsulated drops behave more or less like solid 
particles thus enabling their size distribution measurement 
by a standard size measurement technique. Electronic size 
analysis (Coulter Counter) was employed to measure the 
volume distribution of droplets as it was found to be superior 
to other standard techniques like sieving, sedimentation, 
microscopy and photography. 

5.1.1 Comparison of Electronic Size Analysis with Other Techniques 

Thou^ there are a large number of techniques available 
for particle size measurement, for the size range of droplets 
encountered (50-800p), only a few of the techniques would be 
effective. They are: (i). Sieving (ii) Sedimentation (iii) 
Microscopy (iv) Photography and (v) Electronic Size analysis. 
Figure 5.1 shows the size range of the above methods of 
measurement [60], 




vCoulter principle 
/////// 





un 

L. 
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Madden and McCoy [54] employed wet screening for size 
measurement of the encapsulated drops. They found that nylon 
coating alone is not sufficient for successful sieving as the 
drops tend to become sticky and do not possess enou^ strength 
to withstand the impact. So, they further coated the drops 
with silver thus increasing the strength of the coating. This 
method was tried and was found to he unsuccessful due to 
breakage of encapsulated drops dxiring sieving. 

Sedimentation was tried to obtain the cumulative weight 

fraction for different drop volumes. Water glycerine mixture 

was used as the settling medium in order to increase the 

settling times. The sample was thoroughly mixed by passing 

air from the bottom of the settling column before allowing 

the droplets to settle. This had to be done as it was found 
£ 

to be very difficult to introduce a monolayer of sample at 
the surface of the settling column without agglomeration of 
encapsulated droplets. Two difficulties were encountered. 

(i) Due to high viscosity of the medium, air bubbles continued 
to exist in the column for a long time (ii) Agglomeration of 
encapsulated droplets during settling could not be avoided 
inspite of very low (Oc” per cent) volume fraction of sample. 
Observation of settling droplets throu^ cathetometer revealed 
existence of big ag^omerates which affected the results 
considerably. Mechanical mixing could not be employed as the 
drops were fragile and broke easily. 
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Microscopic counting was also tried for the measurement 
of drop size distribution. Sample was withdrawn from the 
system, diluted and a microscopic slide was prepared. Random 
sites were chosen in the slide and the size and the number 
of droplets were recorded. This was very tedious and time 
consuming procedure involving a lot of personal errors. More- 
over, only the droplet diameters could be measured. Though 
the drops are spherical, a small error in the counts of large 
sized drops gets magnified when the data of number counts 
vs drop diameter is converted to drop volume density as this 
involves evaluation of the third moment of number density with 
respect to drop diameter. 

Until recently, photography was the main techniq.ue by 
which drop size distribution measurements were made. Though 
this technj.lue can measure droplet sizes in situ , it has many 
drawbacks^ It renders only a one dimensional view of the 
droplet population at one particular plane. Possibilities of 
overlapping of droplets and blurring of images of droplets 
not in focus make an accurate drop size measurement very 
difficult. This technique also suffers from the same drawback 
as microscopy in that this also measures only the droplet 
diameter. 

Electronic size analysis measures the volume of droplet 
thus making this technique ideal for the measurement of drop 
volume distribution. Moreover, in situ sampling is possible 
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by placing the probe in the agitated vessel and withdrawing 
the sample soon after encapsulation. This technique is 
very fast and accurate as it is possible to count around 
10,000 droplets for one distribution measurement within 
minutes. 

5.1.2 Principle of Electronic Size Analysis ; 

The number and size of particle in an electrically 
conducting liquid are determined by the application of 0*oulter 
principle [60, 61]. The suspension of particles is forced 
to flow through a small aperture having an immersed electrode 
on each side as shown in Figure 5.2, As each particle passes 
through the aperture, it replaces its own volume of electrolyte 
momentarily changing the resistance value between the electrode 
Since a constant current is maintained between the electrodes 
through a current source, the momentary change in the resistanc 
value between the electrodes produces a voltage pulse of short 
duration having a magnitude proportional to the particle volume' 
Expressing the particle in electrical effect as a 
right cylinder aligned with the aperture axis and shorter 
than the aperture, it can be shown [60, 61] that the change 
in aperture resistance caused by a particle is 

where p = Pesistivity of the medium 
0 

A = Aperture area normal to axis 
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T = Particle voltime 
p = Effective particle resistivity 

a = Area normal to aperture axis of equivalent right 
cylinder for particle as oriented in passage. 
Usually, particle resistivity has been found to be 
many orders of magnitude greater than that of the electrolyte. 
Metal powders and other apparently good conductors behave like 
non-conductors [60]. This is hypothesised as being due to 
oxide surface films and ionic inertia of the Helmholtz electric, 
double layer and associated solvent molecules at the surface 
of each particle [60], In case of drops encapsulated with 
nylon film, one can expect the resistivity of droplets to be 
much higher than the medium, nylon being a good insulator. 

As indicated in equation (5.2,1), deviation from 
linear volumetric response becomes appreciable for nearly 
spherical particles above 50 per cent of aperture diameter 
(a/A - 0,09) [60], This effect is correctable, if need be, ' 
It is markedly reduced for elongated particles, such as fibres,! 
rods and flakes, as the prevailing streamline in the aperture 
causes predominant alignment of such particles with the 
aperture axis. 

The distortion of the electrical field in the aperture 
by a particle will conform to the essential surface of the 
particle rather than follow each crevice and wrinkle, thus 
sensing the ’envelope volume* of the particle. 
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TiiG voltage puLse froin the electrode is amplified, 
scaled and comited to yield the volinne distribution of 
particles. 

5.1.3 Probe Assembly ; 

IPwo glass probes of aperture diameters 2 mm and 0.8 

mm were used. The schematic diagram of the probe assembly 

is shown in Pigure 5.3. Probe consisted of a glass tube 1.5 cm 

diameter fitted with a B-24 female joint, the flattened bottom 

portion of the tube having an aperture of req.uired diameter. 

Platinum electrodes were used. Electrodes were made by ’ 

th 

e.mbedding l/lOOO inch thick platinum foil on a glass plate. i 
The inner electrode ran all along the length of the tube. 

The outer electrode faced the aperture with a gap of 7.5 mm. 

Both the electrodes were kept rigid by connecting the ^ass 
plates carrying the electrodes to the tube. The dimensions 
of the inner and outer electrodes were 11.5 cm x 1 cm and 1 can ; 

1 cm respectively. Here, the distance between the electrodes I 

1 

remained fixed since both the electrodes were embedded on I 

glass plates. This is an improvement on the probe supplied ; 

by Coulter Electronics Inc. for their Coulter Counter Model B 
where both the electrodes were platinum foils hanging into 
the medium. There, there is a definite possibility of ; 

vibration of outer electrode due to the turbilLence of the I 

mediima created by agitation. This changes the distance between' 
the electrodes and gives sp'urious voltage pulses due to the 



|t--*“To vacuum 



Fig. 5-3 -A schematic diagram of the probe assembly! 
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vibration of the outer electrode. Moreover, if the distance 
between the electrodes is not kept fixed then the signal 
strength for a particle of given volume will not be constant 
and hence the results will not be reproducible. 

The probe tube was connected to a vacuum pump supplied 
by Coulter Electronics Inc. through a stopcock. The probe 
assembly was immersed in the suspending electrolytic medium 
and the inner tube was filled with the medium. At the time 
of measurement, suction was applied to the tube by means of 
the vacuum pump and the sample was withdrawn through the small 
aperture by opening the stopcock. The manometer connected 
to the probe gave the pressure drop across the aperture 
during suction, 

5.1.4 Electronic Circuit cry for Measurement ; 

The block diagram of the electronic circuitory is 
shown in Eigure 5.4. The inner platinum electrode of the 
probe was grounded. A positive 60 volts EC was applied between 
the two electrodes through a constant current source and a 
series variable resistor. The current source used was a 
standard current mirror. The current source had provision 
for setting four different currents, namely, (i) 0,5mA (ii) 

1.0 mA (iii) 1,5 mA and (iv) 2 mA, The variable resistor was 
adjusted so as to minimize the voltage drop across the 
transistor in constant current source circuit. The output 
from the probe consisted of a constant D.C, voltage V( =IR,R 




Display 
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being the resistance between the two electrodes) riding 
over which were the voltage pulse signals given by the encapst 
droplets when they pass throng the aperture. The output was 
fed to a passive D.G. filter which filtered the D.C. signal 
thus permitting only the overriding pulse signals to pass 
through. Though shielded cables were used to carry the 
signal and the motors were grounded, there was considerable 
50 cycles pick-up especially from the vacuum pump. Hence it 
was found necessary to use a 50 cycles notch filter to filter 
off any 50 cycles noise picked up by the probe from its 
surroundings. The output from the notch filter consisted only 
of the pulse signals whose amplitude was of the order of millii 
These pulses had to be amplified prior to analysis, BOIL 
pre-amplifier was used for this purpose. Since the amplifier 
had only four discrete gain settings (100, 250, 500 and 1000 ) , 
it was necessary to use an attenuator with attenuations 1, 1/2 
l/5 and l/lO so that a combination of attenuation and gain ; 
settings gave the req.uired amplification. Moreover, large 
voltage pulse signals had to be attenuated before they were 
fed to the amplifier as the output of the amplifier saturated 
at 6 volts. The amplitude of the amplified voltage pulse 
signal was of the order of volts. These signals were very I 

broad with a large rise time (of the order of msecs). They ^ 

had to be shaped before they could be fed to the Multi-Channeli 
Analyser (MCA) for analysis. The pulse shaper converted these; 
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broad voltage piUses into sharp voltage pulses with a very 
small rise time (of the order of microseconds). The amplitude 
of the voltage pulses were preserved during pulse shaping. 

The shaped voltage pulse signals were fed to 128 channel ND-llO 
Multichannel Analyser for analysis. The signals were also 
fed to a Philips Oscilloscope for visual observation. The MCA 
could accept pulses in the range of 0-6.8 volts. The current 
gain and attenuator settings were adjusted so that all the 
amplified and shaped signals fell in this range. The MCA was 
operated in Pulse Height Analysis (PHA) mode. The MCA sensed 
tho pulse amplitude (height) and stored the count of nmber 
of pulses in the appropriate channel number. The MCA was 
calibrated with standard pulses from a precision p^ulse 
generator. The details of the calibration is given in Appendix 
II. Tho calibration gave the relationship between the pulse 
amplitude (height) and the channel number. The MCA gave a 
display of number count versus channel number. It was also 
possible to get a print out of the number count versus channel 
number, 

5.1»5 Size Range of Measurement ; 

Response correction is unnecessary for most parti- 
culate systems if the few largest particles do not exceed 
40-50 per cent of the aperture diameter [60, 61], This is also a 
practical maximum to avoid excessive aperture blockage duie to 
coincidence of large particles. The noise limitations of 
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electronic amplification prevent measurement below 1-2 per cent 
of aperture diameter. Thus, a single aperture size provides a 
diametric measurement range of 20 si, with the corresponding 
volume range upwards of 8000 si [60]. The range of droplet 
sizes encountered in agitated lean liq.uid-liquid dispersions 
was from 50 [j, to SOOp,. Though it is possible to cover the 
whole size range with one probe, in order to increase the 
resolution of volume distribution measurements, it was 
decided to use two probes of aperture diameters 0.8 mm and 
2 mm. The smaller probe was used for the measurement of drop 
population of small sizes - drop population at large stirring 
time and the bigger probe was used for the measurement of 
drop population of large sizesy drop population at small 
stirring time, 

5.1.6 Calibration of the Probes ; 

The amplitude of the voltage ptiLse given by a particle 
of certain volume will depend on many factors such as aperture 
diameter, current, g-ain and a.ttenuator settings and the 
resistance between the two eJ. ectrodes. The relationship 
between the pulse amplitude and particle volume should bo 
known in order to infer the droplet volume distribution from 
the output of MCA, 

A narrow distribution of spherical standard sized 
glass particles were suspended in a conducting medium of haCl 
solution, A sample of the suspended particles was withdrawn 
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through the probe and the spectrum measured by the MCA. The 
average volimae of the narrow sized standard sample was 
measured by microscopic counting. Since the sample is a 
narrow distribution of particles, the average particle volume 
should correspond to the average channel number of the spectrm. 
The same procedure was repeated for different NaCl concentration! 
electronic settings and sample size to obtain the relationship 
between particle volume and channel number for different 
conditions. The details of the calibration is given in 
Appendix II, 

5.1.7 Coinciden ce Eff e ct ; 

Sometimes, more than one droplet may pass through the 
aperture simultaneously. This will give rise to the super- 
imposition of the voJ-tage pulses of these droplets, resulting 
in loss of count, A typical voltage signal from the probe 
when two particles pass through the aperture simultaneously 
is shown in Figure 5.5. As mentioned earlier, the primary 
effect of coincidence is loss of count as some of the 
superimposing voltage pulses will be ignored by the electronic 
circuits. Corrections can be applied to this. These correc- 
tions are known functions of aperture size and particle concen- 
tration and may be checked experimentally by counting a given 
suspension at successive dilutions. As a secondary effect, 
coincidence passages may give rise to voltage pulses of 
amplitude greater than the amplitudes of pulses corresponding 



A typical voltage 

a) One droplet p 

b) Many droplets 
simultaneouslv 
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to either particle thus registering a false particle of 
large volume. This effect will thus vary with the size range 
of the system being analysed. If coincidence loss corrections 
are kept low, this secondary effect will usually be negligible, 
but may be q.uite prc'dominent in narrowly distributed systems 
near the frequency ,peak [60], 

Coincidence corrections were checked experimentally by 
measuring the spectrimi of standard sized spherical glass 
particles of narrow distribution at different concentrations in 
the size range of 50--800p, Since drop volume distribution 
measurements were made at dispersed phase volume fraction of 
0,1 per cent, coincidence effects were verified in the concen- 
tration range of 0,05 to 0.5 per cent. It was found that the 
average channel number corresponding to the standard sample 
did not vary appreciably with concentration indicating that 
the coincidence effects are negligible. 

5.1,8 Bncapsulation Technique; 

Sprow [34] has employed electronic size analysis for 
the measurement of drop size distribution of emulsions by 
withdrawing a sample of emulsion. In case of emulsions, the 
droplet sizes are very small and hence their rate of breakage 
small. Therefore, the drop breakage in the system during 
sampling is very small. The dilution and stabilization of 
the emulsion with glycerine retarded drop coalescence in the 
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system as well as in the aperture tube during sampling. 

However, in case of liquid-liquid dispersions, the droplet 
sizes are fairly large. Therefore, it is necessary to avoid 
breakage and coalescence of droplets during sampling. This 
was achieved by capturing all the droplets in a liquid-liquid 
dispersion in their dynamic state by encapsulating them with 
a thin nylon film. Encapsulation, preserved the identity of 
the droplet by preventing breakage and coalescence during 
sampling. Encapsulation was achieved through an interfacial 
poly condensation reaction. This reaction was almost instanta- 
neous and resulted in the formation of very thin (0,1 - 1 micron 
thick) nylon coating which helped to preserve the identity of 
the droplets at the time of encapsulation thus making a direct 
size measurement possible. Madden and McOoy [543 were the 
first to illustrate the use of this technique for the measure- 
ment of drop size distribution. Encapsulation was employed by 
Mylnek and Resnick [39] for the measurement of drop size 
distribution in liquid -liquid dispersions. 

The dispersed phases used for this study were CCl^-i- 
Octane mixture (50—50 per cent), Anisole— CCl^ mixture (80—20 
per cent) and chlorobenzene. The reactants for the interfacial 
poly condensation reaction were terephthalic acid chloride and 
piperazine [39], Terephthalic acid chloride was present in 
the dispersed phase at a concentration of 1 per cent by weight. 
At the time of encapsulation, mixture of piperazine and HaOH 
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solutions "was added to the continuous phase for the roaction 
to occur. The reaction was found to he instantaneous resulting 
in a coating of Hylon over all the droplets. The nylon 

film, being very thin ( 0,1 - l.Op), did not affect the 

drop volume measurement. The details of the reaction and the 
composition of the reactants are given in Appendix III, 

5.2 EXPERIiffiUTAL SET-UP; 

The details of the mixing vessel assembly are given in 
Figure 5.6, The mixing vessel was a 5 litre jacketted glass 
vessel, 14 cms in diameter and 20 cms high. It was filled 
with the liq.uid-liq.irLd mixture to a height of 16,25 cms. The 
vessel did not have any sharp corners at the bottom to avoid 
dead zones. The l/D ratio was very nearly one as per the 
recomiiiendations of Shinnar and Church [28], Outer jacket was 
provided for water circulation from a constant temperature 
water bath in order to maintain a constant temperature. All 
experiments were conducted at 30 + O.l^C. 

The mixing vessel was placed on a bottom aluminium 
plate fixed to the work bench. Rubber gasket was provided 
between the vessel flange and the stainless steel baffle 
assembly. The baffle assembly consisted of a top stainless 
steel plate and four baffles. The baffles were 19 cms long, 
1,4 cm wide and 0.625 cms thick. The baffle assembly was 
connected to the bottom plate through three brass rods and 

was kept rigidly in position by tightening the brass flynuts. 

A teflon bush was connected to the top stainless steel plate. 
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Fig. 5-6-Plan and elevation of mixing vessel assembly. 
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This acted as a guide to the agitator shaft and avoided 
wobbling of the agitator. The number and dimensions of the 
baffles were fixed based on 'Standard baffling'. Fernandes 
and Sharma [62] used four vertical baffles, each 1/lOth of 
the vessel diameter, mounted against the vessel wall at right 
angles to it, spaced at 90° inteivals around the vessel. Mack 
and Kroll [63] and Rushton, Costich and Evere1:t [15] have 
shown that four baffles are satisfactory for 'standard bafflin 
According to Mack and Kroll [63], baffles slightly away from 
the wall would give more power consumption and hence better 
mixing than when the baffles are attached to the wall. On the 
basis of these considerations, four baffles each l/lOth of 
the vessel diameter were used. The baffles were equally 
spaced at 90°, right angles to the vessel and 1/4' inch away 
from the wall. 

Two ordinary 6-blade paddle impellers of 2 inch and 3 
inch diameters were used. Most of the experiments were conduo 

with 3 inch diameter impeller to ensure good mixing. The 

th "th 'i 

impellers were 3/8 inch wide and 1/16 inch thick. The impellej 

was placed at 4 cm from the bottom of the vessel. 

The agitator was driven by a 1/4 H.P, induction motor 
with a variable speed drive. The variable speed drive 
consisted of a combination of step pulleys and gears as the 
speed reduction ratio required was high. Experiments were , 
conducted at three speeds - 300 RPM, 420 RPM and 480 RM. ! 
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It -was not posylble to con.dinct the experiments at higher 
speeds as the drop volLiine distribution changed rapidly due 
to rapid drop breakage, 

5.3 EXPERIMENTAL PROCEDURE; 

The mixing vessel, the stainless steel baffle assembly 
and the agitator ■were cleaned with chromic acid and rinsed, 
thoroughly with distilled water. The mixing vessel was 
placed on the bottom aluminium plate. The baffle assembly 
was centred in the vessel and rigidly fixed by tightening the 
brass flynuts. The agitator was connected to the variable 
speed drive. The height of the impeller above the vessel 
bottom was checked. The probe assembly was lowered into the 
vessel through a hole provided in the top stainless steel 
plate so that the apertu.re of the probe was placed between two 
baffles and at half the depth of the vessel. The vessel was 
then charged with 2500 cc of distilled water presaturated 
with the dispersed phase. Presaturation was done to avoid 
mass transfer of the dispersed phase into continuous phase 
during agitation, ¥ater was circulated in the outer jacket of 
the vessel from a thermostat. The system was allowed to 
eq.uilibrate for abou.t half an hour, after which it was found 
that the continuous phase temperature reached the set point 
in the thermostat. The continuous phase temperature was 
measured with a thermistor probe using a tele thermometer . All 
experimental measurements were made at 30 + 0,1°C. 
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One pe^:* cent by weight of terephthalic acid chloride 
was dissolved in the dispersed phase. 2,5 cc of this was added 
to the continuous phase maintaining a dispersed phase volime 
fraction of 0.1 per cent. The induction motor was started 
and the time of agitation was measured with a stopwatch. The 
speed of agitation was checked for every observation with 
General Radio Strobotac. At the time of measurement, 50 cc 
of mixture of piperazine and NaOH was added to the continuous 
phase and stirring continued for 10 seconds after which the 
speed of agitation was reduced. The dispersion was agitated 
at lower ^eed for a few minutes for the interfacial poly- 
condensation reaction to be completed in order to increase 
the strength of nylon coating. Because of the presence of 
RaOH and formation of NaCl in poly condensation reaction, the 
continuous phase became a conducting medium, Vacuum was 
applied to the probe assembly throu^ the vacuum pump. The 
stopcock in the outlet of the probe tube was then slowly 
opened to fill the probe with the dispersion. The resistance 
between the electrodes was then measured with a multimeter. 

The current, gain and attenuator settings were adjusted. 
Accordingly, the resistance of the series variable resitor was 
also adjusted, ' Electrophoresis power supply was put on to 
apply 60 volts DC between the electrodes through the constant 
current source. The output from the amplifier was observed 
in the Philips oscilloscope. How, a sample of the dispersion 
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■was 'with.drawn through the aperture hy slowly opening the 
stopcock at the outlet of the probe. The pressure drop 
across the aperture "was checked with the mercury manometer 
connected to the probe. The pressure drop was adjusted by 
adjusting the opening of the stopcock. Same pressure drop 
was maintained during sampling for all observations. This 
was done to ensure that enough vacu'um was applied to withdraw 
large size droplets. The amplified voltage pulse signals 
could be observed in the oscilloscope. The voltage range of 
the pulses was noted. If the pulses were not in the proper 
range (either too small or too big) then the current, 
attenuator and gain settings were adjusted accordingly so as 
to get the pulses in the proper range. The output from the 
pulse shaper was observed in the oscilloscope. If found 
satisfactory, it was fed to the MCA and the pulses were logged 
in the MCA in PHA mode. Usually, a sample of 200 cc was 
withdrawn at a time. Then the spectrxmn was observed in the 
visual display of the MCA, If the total n'umber count of the 
droplets was not sufficient, then the withdrawn sample was 
recycled to the vessel and another sample of 200 cc was with- 
drawn, pulses logged in MCA. Around 10,000 droplets were 
counted to obtain a vol'ume distribution. The printout of the 
spectr-um was then taken. The same procedure was repeated to get 
another reading of the volume distribution. This was done to 
check the consistency of the measured distribution. The voltage 
drop between the electrodes was then measured with the 
multimeter. 



CHAPTER 6 


RESULTS MB DISCUSSIOF 

This chapter presents the experimental results of 
transient and eq,uilibrium drop volume distributions in batch 
agitated lean liquid -liquid dispersions. The experimental 
results are also compared with the models for transitional 
breakage probability and equilibrium drop size distribution. 

6.1 SYSTEMS MB EXRERIl^IENTAL GONDITIOHS; 

In all the experiments the continuous phase was 
distilled water presaturated with the dispersed phase liquid. 

Drop volume distribution measurements were made for three 
different dispersed phases and three agitator speeds. Three 
inch diameter paddle impeller was used in most of the experi- 
ments, Measurements were also made with two inch diameter paddle 
impeller in order to study the effect of impeller diameter. 

In all the experiments, the density and viscosity of both the 
continuous and dispersed phases were more or less equal. In 
fact, this dictated the choice of the dispersed phases. All 
the experiments were conducted at 30 ± 0.1°C and dispersed phase 
volume fraction 0 was maintained at 0.1 per cent in most of 
the experiments. Measurements were also made at three more 
dispersed phase volume' fractions, namely 0.25 per cent, 0.375 
per cent and 0.5 per cent, to study the effect of 0 on drop 
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voliime distribution. The different systems and the experi- 
mental conditions are given in Table 6.1. The interfacial 
tension of the three systems were measured at 30 + 0.1°C by 
drop weight method. The physical properties of the dispersed 
phases are given in Table 6.2. At least 10,000 droplets 
were counted for each drop volume distribution measurement. 

In all the measurements, the sample was withdrawn from the 
same position in the vessel. The sampling probe was positioned 
to keep the aperture exactly at half the liquid depth midway 
between two baffles, Uo attempt was made to study the effect 
of position on drop volume distribution as it was found that ^ 
the variation of drop size distribution with position is 
negligible in systems where the buoyancy effects are negligible 
[ 39 ]. 

6.2 REPRODUCIBILITY OE MEASUREMENT: 

Since drop breakage and coalescence are random in agi- 
tated liquid -liquid dispersions, it is necessary to count 
large number of droplets to get a meaningful distribution data. 
Therefore, at least 10,000 droplets were counted for each drop 
volume distribution measurement. A typical plot of number 
count of droplets versus channel number is shown in Figure 6,1, 
In each experiment , at least two samples were withdrawn from 
the drop population in order to ascertain that the samples 
were representative of the population, Kolmogorov-Smirnov 
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TABLE 6.1; SYSTEMS M3 EXPERIMBHTAL COEDITIORS 


System 

Agitated Speed 
RPM 

Impeller Dia, 
inches 

Wat er-CCl.+i -Octane (50-50 
^ per cent ) 

480 

3 and 2 

¥ater-CCl,+i -Octane (50-50 
^ per cent) 

420 

3 

Wat er-CCl.+i -Octane (50-50 

per cent) 

300 

3 

Water-Anisole+CCl^ (80-20 

per cent) 

300 

3 

Water-Chlorobenzene 

300 

3 

TABLE 6.2: PHYSICAL PROPERTIES OP THE DISPERSED PHASE 

Dispersed Phase 

Density , gms/ cc 

Interfacial 
tension, dynes/ ' 

OCl^+i -Octane (50-50 per cent) 

1.149 

46 

Anisole+CCl^ (80-20 per cent) 

1.113 

29 

Chlorob enzene 

1.101 

37.7 




Number count / channel 



Fig. 6-1 -Plot of number count of droplets versus 
channel number. 
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two sample test was employed to verify that the two samples 
represented the same population, A typical comparison of the 
cumulative number per cent of two samples from the same 
population is shown in Figure 6.2, The result of the 

Kolmogorov-Smirnov two sample test is given in Table 6.3, 

Kolmogorov—Smirnov two sample test was found to be very 

stringent as the total number count of droplets for each 

measurement was very large. The test showed that the 
probability that the two samples are from the same population 
is finite though very small. Experiments were also repeated 
under identical conditions to verify the reproducibility of 
the experiment. The reproducibility was verified by means 
of Kolmogorov~Smirnov two sample test, A typical comparison 
of the cumulative number per cent of two experimental obser- 
vations is made in Figure 6.3. The result of Kolmogorov- 
Smirnov two sample test is given in Table 6,3. Thou^ the 
experimental cumulative number per cents are not very different, 
Kolmogorov-Smirnov test, being very stringent, gave a very 
small probability that the two samples are from the same 
population. 

6.3 EFFECT OF DISPERSED PHASE FRACTIONS 

At low dispersed phase fraction coalescence should be 
negligible due to negligible interaction between droplets. 
Therefore, one would expect the effect of dispersed phase 
fraction 0 on drop volume distribution for lean agitated 
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Fig. 6-2 -Comparison of cumulative number percents of 
two samples. 
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TABLE 6.3 ; BE SUITS OF KOMO&OROY-SMIRITOV TWO SAMPLE TEST 


EEI«I-300 ¥ater-C01.+i -Octane (50-50 

^ per cent) 


Probability that both 
the samples are from th 
same population (per cej 


Time-15 min. 3 in paddle Temp,, 30 0 
Sample 1 
Sample 2 

BPM-300 ¥ater-Anisole+CCl^ (80-20 

per cent) 

Time-1 min 3 iu, paddle Temp,, 30^0 


1.325 


Experiment 1 


0.1138 


Experiment 2 


CumutativG number percent 
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Drop volume (cc) 


Fig. 6-3 - Comparison of cumulative number percents 
of repeated experiments. 




115 


liq.uid-llq.uid dispersions to le negligible. Drop volume 
distributions were measured at three more dispersed phase 
volume fractions 0.25 pe^ cent, 0,375 cent and 0,5 per 
cent at the agitator speed of 480 RM for water - CCl^ + 
i-octane (50-50 per cent) system, Pigure 6,4 compares the 
cumulative number per cents at different dispersed phase 
volume fractions, Kolmogorov-Smimov two sample test showed 
that the distributions represent different popiilations. Prom 
Pigure 6.4, it can be seen that the cumulative number per 
cent shifts to larger drop volume as the dispersed phase 
volume fraction increases, though the shift is very small. 
This shows that there is a definite but negligible effect of 
dispersed phase volume fraction on drop volume distribution. 
Hence, in all the experiments it can bo assumed that coales- 
cence is negligible. 

6.4 TRANSIENTS OP DROP YOllME DISTRIBUTIONi 

Transients of drop volume distribution were measured 
for water-CCl ^-♦'i-octane (50-50 per cent) system at three 
different agitator speeds of 480 RPM, 420 RPM and 300 RPM 
with three inch diameter paddle impeller and at the agitator 
speed of 480 RPM with two inch diameter paddle impeller. 
Transients of drop volume distribution were also measured at 
the agitator speed of 300 RPM with three inch diameter paddle 
impeller for water-Anisole + CCl^ (80-20 per cent) and water- 
chlorobenzene systems. Pigures 6.5 throu^ 6.10 show the 
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Fig. 6-5 : Effect of stirring time on drop volume distribution 
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effect of stirring time on drop volume distribution for the 
above measurements. 

6,4.1 Verification of * Power Law* Approximation of Transitiona! 

Breakage Probability : 

If the ‘power law* approximation Kv^ for the transi- 
tional breakage probability T (v) is valid, it should be 
possible to collapse the drop volume distributions at differen' 
times into a single curve under the similarity transformation , 
v^t for sufficiently large time [4]. Moreover, log-log plots, 
of droplet volume v versus time t for different fixed values 
of cumulative volume per cents should yield a set of parallel 
straight lines of slope - ~ [4]. Piguresb.ll through 6.16 
show plots of droplet volume v versus time t for different 
fixed values of cumulative volume per cents. In Figures 6.11 
through 6.16, it was found that the data could be represented 
by two sets of parallel straight lines of different slopes 
i.e. the transitional breakage probability r(v) could be 
approximated by a ‘power law’ of varying exponents. The sets 
of data in Figures 6.11 throu^ 6.16 could be divided into 
two regions of different exponents. The region corresponding 
to large droplet volumes or small times will be referred to 
hereafter as region I and the other as region II, Earlier, 
it was shown from Figure 3,2 that the transitional breakage 
probability predicted by the model could be approximated by 
a series of ‘power law' models with different exponents for 
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Fig. 6-13 - Plot of stirring time vs. droplet volume for a fixed cumulative 
volume percent. 




Fig. 6-14 -Plot of stirring time vs. droplet volume tor a fixed cumulative 
volume percent. 




Fig. 6-15 ; Plot of stirring time vs. droplet volume for a fixed cumulative 





Fig. 6-16 ; Plot of stirring time vs, droplet volume tor o fixed 
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different regions of droplet ■volume 'with, the value of the 
exponent increasing ■with the decrease in droplet volume, 
figure 3,2 gives the value of the exponents for different 
regions of droplet volumes. Thou^ all the regions could not 
be identified in one set of data, in all the cases two distinct 
regions of different value of the exponent co'uld be detected. 
Moreover, the value of the exponent n obtained from the 
slope of the set of parallel straigh't lines was found to 
increase at large times, or eq.uivalentlyj for small droplet 
volumes. Tables 6,4 through 6,9 present the values of the 
exponent *n’ obtained through least square fits for different 
cum'ulative vol'ume per cents for both the regions. Sets 
of more or less parallel straight lines were obtained in all 
the cases especially in the cum'ulative volume per cent range 
of 40-70 per cent. However, there is deviation for small 
and large values of cum'ulative volume per cents. This cculd 
be attributed to the experimental error. All the drop volume 
distributions sho'uld exhibit zero slope at 100 cumulative 
volume per cent. But, some drop vol'ume distributions do not 
exhibit zero slope at 100 cumulative volume per cent as can be 
seen from I’igures6,5 through 6,10. This was due to loss of 
co'unt of some large voltage pulses corresponding to large 
vol'ume droplets. With the available combinations of discrete 
gain and current settings, it was sometimes not possible to 
choose the required amplification without losing either the 



131 


TABLE 6.4 s V AL’JBS OF THE » POVfER LA¥» EJCPONEBT FOR LIFEEREM} 

CUMULATIVE' VOLUI4E PER CENTS 


RPM - 480 ¥ater-CCl^ 

+i -Octane (50-50 per cent) 

Temp.30°C 

Impeller; 3 in, paddle 

Exponent ’n' 


C-um.Vol. Per cent 

Region I 

Region II 

20 

1.276 

1.766 

50 

1.186 

1.783 

40 

1.006 

1.920 

50 

1.096 

2.226 

60 

1.058 

2.048 

70 

1.066 

2.189 

80 

1.039 

2.135 

TABLE 6. 5 5 VALUES OP 

THE ’POV/ER LAW EZPOFENT EOR DIPEBRENT 


CMULATIVE VOLUME PER CENTS 

RPM 420 Water-CCl^+i 

-Octane (50-50 per cent) 

Temp. 30° C 

Impeller. 3 in. paddle 

Exponent ’n' 


Oum.Vol, per cent 

Region I 

Region II 

20 

1.314 

1.577 

30 

1.191 

1.569 

40 

1.091 

1.580 

50 

1.044 

1.559 

60 

1.004 

1,568 

70 

0.978 

1.581 

80 

0.944 

1.592 
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TABLE 

6.6; 

VilUES OE 'POWER LAW EXPOlffilTT FOR : 

DIFFERENT 

Rm - 300 


OIMDLATIYE YOLJME PER CEElS 

Wat er~CCl^+i -Octane (50-50 per cent) 

Temp . 30°< 

Impeller 

- 3 

in. paddle 


Cxun.Tol, 

per 

Exponent 

cent Region I 

'n' 

Region Ii 

20 


1.213 

1.288 

30 


1.160 

1.274 

40 


1.096 

1.362 

50 


1.107 

1.365 

60 


1.081 

1.486 

70 


1.068 

1.483 

80 


1.050 

1.556 

TABLE 

. 6.7 

VALUES OE 'POWER LAW EjCPOMT FOR 

DIFPEPIENT 

RPM - 300 


CIMULATIVE VOLUME PER CENTS 

¥ater-Anisole+CCl^ (80-20 per cent) 

Temp , 30° 

Impeller 

- 3 

in. paddle 


Oum.Vol, 

per 

Exponent ' 

cent Region 1 

n' 

Region H 

20 


1.113 

1.894 

30 


1.060 

1.888 

40 


1.039 

1.789 

50 


1.030 

1.737 

60 


1.028 

1.667 

70 


1.036 

1.578 

80 


1.046 

1.539 
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TABLE 6.8: 

VALLES OE THE 'POWER LAW 

' EIPONEHT 

POR DIPPEREHT 

EPM 300 

CLMULATIVE VOLUIffi 

Water-Chlorobenzene 

PER CEHTS 

Temp . ; 30° C 

CuEL, Vol. per 

Impeller; 3 in. paddle 

cent Region 

Exponent ',n' 

I Region II 

20 

1.900 


2,190 

30 

1.674 


1.946 

40 

1.551 


2.003 

50 

1.465 


1.893 

60 

1.452 


1.877 

70 

1.455 


1.615 

80 

1.433 


1.683 

TABLE 6.9; 

VALUES OP THE 'POWER LAW 

' EXPOHENT 

POR DIPFEREIT 

Cuin. Vol, per 

CUMULATIVE VOLUME 

cent Region 

PER CEHTS 

Exponent 

I 

'n' 

Region II 

20 

1.370 


2.267 

30 

1.285 


2.044 

40 

1.253 


1.811 

50 

1.224 


1.588 

60 

1.199 


1.467 

70 

1.169 


1.359 

80 

1.148 


1.332 
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sensitivity or some piiLses of large amplitude. 

Excepting for water— chlorobenzene system, the value 
of the exponent n for region I was around 1,0, Eor water- 
chlorobenzene system, the value of n for region I was around 
1,5 • Ihe value of n for region II for water-CCl^+i -Octane 
(50-50 per cent) system at the agitator ^eed of 480 RPM 
was around 2.0. Eor water-chlorobenzene and water-ihiisole + CC 
(80—20 per cent) systems at the agitator speed of 300 RPM 
the value of n for region II was around 1,8. In all the other 
cases, the value of n for region II was arotind 1,5. Hence, 
three distinct regions of e^iponents 1.0, 1.5 and 2.0 could be j 
identified. The intermediate value of 1.8 obtained for the . 
exponent n in two cases is believed to be due to overlapping of 
two regions. It would be interesting to note that in case of 
low energy input i.e, at the agitator speeds of 300 RPM and 
420 RPM and at the agitator speed of 480 RPM with two inch 
diameter paddle impeller, the region with the exponent of 1,5 
could be identified. The value of the exponents obtained 
from the data agree fairly well with the values of 1,08, 1.42, 
and 2,12 obtained from the approximation of the predicted 
transitional breakage probability by a series of 'power law' 
models of varying exponents. The regions with exponent 0.574 
and 3.0 could not be identified in any of the experimental 
data. This is understandable as drop volume distribution 
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measurements were not made either at ■very small times or 
very near e(iuilibri-um. 

6.4.2 Verification of the Model for Transitional Breakage 
Probability ; 

As the transitional breakage probability r(v) predicte 
by the model can be approximated by a series of 'power law* 
models Ev'^ with different values of the exponent n for 
different droplet volume regions and for each of these regions 
Ev .t is the similarity transformation [4] which would collaps, 
all the drop volume distributions in that region into a single 
curve, one would expect that a transformation r(v)*t would 
collapse all the drop volume distributions at all times into 
a single curve irrespective of their regions if the model 
for the transitional breakage probability ■ is indeed valid. 
However, the computation, of transitional breakage probability 
as predicted by the model, 

r (v) = X erfo [ { 

^-5/18/2 - 1/3] 

yielded a very large value of the maximum stable drop volume 
Vg, In other words, for most of the drop volume range 
encountered in the experiments, the computed transitional 
breakage probability turned out to be zero. For the compu- 
tation of transitional breakage probability, power input per 
unit mass e was calculated from the empirical relation. 
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e = 5.1 IT^L^ (6.4.1) 

based on the experimental correlation [15 » 16] for 2 blade 
paddle impeller e = 1,7 under the assimiption that the 

power input is proportional to the number of blades, A factor 
was therefore introduced into the argument of the complimentarj 
error function in order to match the variation of F (v) with 
droplet volume v in the drop volume range encountered with 
the experimentally observed variation, This had to be done 
by trial and error, A factor of 0,2 was found to represent 
the variation of transitional breakage probability r (v) with 
droplet volume v satisfactorily, 

Prom the observations of Cutter [18], we find that 
about 20 per cent of the energy is dissipated in the impeller 
itself, about 50 per cent in the impeller stream, about 30 
per cent in the other parts of the stirred vessel. Measure- 
ments of s/e in the impeller stream showed' an almost expo- 
nential variation from 70 near the impeller tip to about 
3.5 near the wall. However, Gimkel and Weber [24] calc\xlated 
the energy dissipation rate e in the impeller stream throu^ 
energy balance and showed that most of the energy is dissi- 
pated in the bulk of the vessel, Ihey also claimed that the 
observations of Cutter were due to experimental errors. If 
such important inhomo geni ties do exist in an agitated vessel 
as claimed by Cutter [18] , it is felt that some statistical 
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mean e should be used instead of e in the probability 
distribution of relative velocity. Therefore, a statistical 
mean e should be used instead of e in the expression for ^(v) . 
It may be reasonable to use e averaged over the residence 
time of the drops in different regions of the stirred vessel 
i.e, in effect, the mean s may be the average energy 
dissipation rate a drop 'experiences’ in the stirred vessel. 
This will depend on the flow pattern in the agitated vessel. 

If the inhomo geni ties in the stirred vessel are considerable, 
one would expect most of the breakage to occur in the impeller 
stream. In such a case, the dissipation rate per unit mass a 
droplet 'experiences' will be much larger than e. This may be 
the reason for the prediction of large maximum stable drop 
volume Vg from the model for I (v) using e, thus necessitating 
the introduction of the factor 0,2, It may also be possible 
that the energy barrier for breakage of a droplet of volimie , 
V, i.e, the minimum increase in the inter facial energy of a 
droplet of volume v for fragmentation, is much less than the 
increase in the interfa-cial energy for binary equal breakage. 
Hence , the experimental drop volume distributions were 
transformed using the transformation erfc[{ ^ 

(^)l/9 ^ e Since the factor 0,2 introduced 

in the arguement of the complimentary error function was 
obliained by matching the variation of computed r(v) with 
the experimental variation for three inch diameter paddle 
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impeller and tlie value of th.e factor for two inch diameter 
paddle impeller is believed to be different due to the 
difference in the flow pattern in the stirred vessel, the 
transformation was performed only for the experimental drop 
volume distributions for three inch diameter paddle impeller, 
Figure 6,17 shows the plot of transformed drop volume 
distributions. Most of the transformed drop volume distri- 
butions seem to fall on a single cuive though there is 
considerable scatter especially in case of drop volirme distri' 
butions at the agitator speed of 300 REM. 

The transients of the drop volume distribution were 
found to agree fairly well with the ’power law’ model for 
transitional breakage probability with varying exponents. 

The values of the exponents 1,0, 1.5 and 2,0 obtained for 
the three regions identified experimentally also agree fairly 
well with the values of 1.08, 1,42 and 2.12 obtained from 
the approximation of the predicted transitional breakage 
probability by a series of 'power law’ models. The trans- 
formation 

erfo [ { 2{2W3_i) 2 (|)l/9 ^-SAS/iq i 

seems to collapse all the drop volume distributions into a 
single curve. These findings validate the model for transi- 
tional breakage probability. 
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Fig 6 17 - Plot of cumulative volume percent versus 
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6 *4.5 Estimation of Daughter Droplet Distribution ; 

The drop ’volume distributions in a region “were trans- 
formed using a modified similarity transformation 

_ "V 

V and n being an average droplet vol'ume and the exponent of 
’power la'w’ approximation of r (v) for that region respectively 
Figures 6,18 through 6,29 show the plot of the modified drop 
vol’ume density f { versus (^)*t4/^ for regions I 

V V 

and II, It can be seen from Figures 6,18 through 6.29 that the 

modified drop vol’ume densities do collapse into a single curve. 

There is considerable scatter at very lo’W’ droplet volumes 

especially at lo’w energy input. This can be attributed to 

the experimental errors. The solid c’urves in Figures 6,18 

through 6,29 are the modified drop volume densities obtained 

through loast-s<luare fit of the data using polynomial expansions'. 

2 

The error s’um, of squares X for that fit are given in Figures 

p 

6,18 through 6.29« X was found to be much larger for region I 
indicating more scatter. This was due to more experimental 
error in the measurement of drop volume distributions at small 
stirring times. 

Since the similarity transformation collapses 

V 

all the drop volume densities in a region into a single curve, 
the assumption of 'similar breakage' made for the derivation 
of the similarity transformation [53] seems to be valid. 
Therefore, the moments of the function Kg(x) could be 
computed from the relation [4]» 
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Fig. 6-20- Plot of modified drop volume density 
versus similarity variable (v/v)t^^*^. 
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Fig. 6- 22 - Plot of modified drop volume density f|(v/v)t 
• similarity variable (v/v)!^^*^. 


4 


m 


RPM-300 T(zmp.-30T 
Impeller - 3 in paddle 
Continuous phase- Water 
Dispersed phase- CCI 4 + i -octane 

(50-50°/o) 

Region II 
v=0' 2118 X lO^cc 
x2^0.2611 xlO'^ 

o 1 5 m i n 
• 20 min 

’ ^ 25min 


V(ftcP 


• 


Fig 6-23-Plot of modified drop volume density f|(v/v)t 
versus similarity variable (v/v) t"'^'^. 
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Fig. 6-24- Plot of modified drop volume density f|(v/v)t 
versus similarity variable (v/ v) 
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Fig, 6-26 -Plot of modified drop volume density f|( v/ v) 
versus similarity variable (v/v) t^^^ . 
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Fig. 6-29 - Plot of modified drop volume density f|(v/v)t '^'^ 
versus similarity variable (v/ v)t^^'^. 





153 


^r+n ^ 

where g( 2 c) is the normalised daughter droplet distribution 

function for ’similar breakage' and 

X = J d »? Z being the similarity 

0 

transformation v t^^^. 

The moments { X^ } were computed from the moments } 
of the modified drop volume density f{ (“)*t^^^} from the 


relation, 


V 




The modified drop volimae density obtained through least square 


fit of the data was used for the calculation of the moments{Xp}. 

The function Kg(x) was then estimated from its moments 
expressing the function in terms of orthogonal poly- 
nomials generated from a suitable weighting function. The co- 
efficients of orthogonal polynomial expansion were estimated 
from the moments { } • Three weighting functions were used 

for the estimation of the function Kg(x) in order to ascertain 
the consistency of estimation. They were;(i) 1 (ii) x and (iii) 
3x — 2x'^ , an approximation of gaussian distribution. In case of 
the first two weighting function^, the convergence was found to 
be poor even after inclusion of four terms in the expansion^ 


Inclusion of more number of tenns in the orthogonal expansion 
was desisted as the evaluation of higher coefficients would 
require the knowledge of hi^er moments of the modified drop 
volume density f and the errors involved in the 

V 

estimation of higher moments X™ are considerable. Hence, 
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•the fimction iCg(x) could no'b be estimated with the first two 
weighting f'unctions, Figure 6.30 shows the estimated function 
Kg(x) through orthogonal polynomial expansion with the weighting 
function 3x^-2x^ for 2 <3 and 4 terms in the expansion. From 
Figure 6,30 it can be seen that the convergence is satisfactory. 
Four terms were included in the orthogonal polynomial expansion 
for the estimation of the function Eg(x). The rate constant Z 
for the ‘power law’ approximation Kv^ of the transitional 
breakage probability was then estimated from the function Kg(x) 
by making use of the constraint g(l) = 1. Table 6.10 gives 
the values of the computed rate constants. Table 6,11 gives 
the computed moments { normalised cimiulative 

daughter drop distribution fiznction g(x) . It can be seen from 
Table 6,11 that the moments of g(x) decrease progressively with ; 
the increase in the order of the moments. This is understandable 
as g(x) is non-zero only in the unit interval. In case of 
water- Inisole + CCl^(80-«20 per cent) and water-chlorobenzene 
systems the computed moments of g(x) for region II were found 
to be high. In all the other cases the value of the monienbs 
compare fairly well with a variation of around 20 per cent. 

From the computed values of the moments of g(x) it seems that 

I 

the normalised c-umulative daughter droplet distribution 
function is more or less same for different systems and experi- | 
mental conditions. This is further indicated by the plots of : 
estimated cumulative daughter droplet distribution f-unction g(x) | 
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T ABLE 6.10 ; VALUES OE RATE COKSTAETS 


System 

Agitated 

SpeedjEPM 

Impeller 
dia, , in. 

Region I 

Region II 

Water-CCl.' 
(50-50 per 
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cent) 
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and the density g’ (x) shown in Figures 6,51 through 6.42, 
from which it is evident that the nature of the dau^ter 
droplet distribution function is the same in all the cases. 

The above results indicate that the mechanism of drop breakage 
is the same for different experimental conditions and systems. 
Since g(x) is the emulative daughter droplet distri- 
bution function, its zeroth moment would determine the average 
daughter droplet volume and its first moment would yield the 
standard deviation about the mean. Table 6.12 gives the 
average non-dimensionalized daughter droplet volume and the 
standard deviation for all sets of experimental data. The 
values of average daughter droplet volume obtained for 
different sets of experimental data compare fairly well 
whereas there is appreciable variation in the values of 
standard deviation. The standard deviation was found to be 
fairly high which indicates that the daughter droplet density 
is quite broad. Table 6,13 gives the mean and standard 
deviation of the dau^ter droplet distribution g(x) estimated 
through orthonormal polynomial expansion. From Tables 6,12 
and 6,13, we find that the mean daughter droplet volume 
computed from the moments of g(2:) does not compare well XArith 
the mean daughter droplet volume of the estimated g(x) 
whereas the standard deviation in both the cases compare 
fairly well excepting for water-chlorobenzene and water-Anisole 
+CCl^(80-*20 per cent) systems for region II, The mean daughter 



Cumulative volume fraction g(x) 
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Normalised daughter droplet • volume x 


Fig. 6-31 -Plots of estimated daughter droplet distribution 
and density . 


Daughter droplet density g'(x) 
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Fig. 6-32- Plots of estimated daugher droplet distribution 
and density. 


Daughter droplet density g'(x) 
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Fig. 6-33 - Plots of estimated daughter droplet distribution 
and density. 


Daughter droplet density g'(x) 



Cumulative volume fraction g(x) 



Cumulative volume tracti 



Fig. 6-35 - Plots of Gstimatcd 
and density. 


daughter droplet distribution 


Daughter droplet density g'(x) 


Cumulative volume fraction g(x) 
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Normalised daughter droplet volume x 


Fig. 6-36 -Plots of estimated daughter droplet distribution 
and densitv . 


Daughter droplet density g'(x) 
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Fig. 6-37 - Plots, of estimated daughter droplet distribution 
and density . 
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Fig. 6-38 - Plots of estimated daughter droplet distribution 
and density. 


Daughter droplet density g'(x) 



Cumulative volume fraction g (x) 


167 



Fig. 6-39 - Plots of estimated daughter droplet distribution 
and density . 


Daughter droplet density g'(x) 
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Daughter droplet density g^(x) 
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Fig. 6-41 - Plots of estimated daughter droplet distri bution 
and density. 


Daughter droplet density g'(x) 
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Fig. 6-42-Plots of estimated daughter droplet distribution 
and density. 


Daughter droplet density g'(x) 
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droplet volume of the estimated gix) is found to be lower than 
the moan daughter droplet volume computed from the moments 
of g( 2 :) , Moreover, the estimated daughter droplet density 
became negative for large daughter droplet volumes when the 
standard deviation computed from the moments of gCx) is high. 

The estimated daughter droplet distribution and density were 
found to be particularly bad in case of water chlorobenzene 
system for region II. It can also be observed that the value 
of estimated daughter droplet density at normalised daughter 
droplet volume of 1,0 was found to be a large positive value 
in all cases where the standard deviation computed from the 
moments of g(x) was small. In view of these, it may be 
necessary to include more number of terms in the orthogonal 
polynomial expansion for accurate estimation of function Kg(x). 
As stated earlier, such an attempt was desisted because of 
large errors involved in the calculation of higher moments. 
Therefore, the estimated daughter droplet distributions are 
only qualitative estimations. 

The following inferences could therefore be made from 
the experimental results: 

(i) The assumption of 'similar breakage' is valid. 

(ii) The breakage mechanism seems to be the same for a 
all the systems and experimental conditions, 

(iii) The noimalised daughter droplet density is 
more or less same for all the cases. 
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(iv) Tke daughter droplet density is broad, 

6.4.4 Variation of the Rate Constant K with e and a t 

The functional dependence of the rate constant K in 
the ’power law’ approximation of the transitional breakage 
probability on the physical properties of the system and the 
experimental conditions is given by, 

n 18n 12n ~ 9n 

Zav”alT5 o 5 

a 

Figure 6,43 shows a log~log plot of the rate constant K 
versus the agitator speed for water - OCl^+i -Octane (50-50 
per cent) system for region I. Since the value of the 
exponent n for region I was around 1,0, the functional 
dependence of the rate constant K on the agitator speed N is 
according to the model, 

Z a 

From Figure 6*42 it can be seen that the computed values 
of Z fall on a straight line of slope 3.38, Hence, the 
functional dependence of the computed values of Z on the 
agitator speed Z for region I agrees fairly well with the 
model. For region II, the value of the exponent n for 
water— CCi^ +i— ociiane (50—50 per cent) system at the agitator 
speeds of 300 RPM and 420 RPM were around 1,5 and at the 
agitator speed of 480 EPM the exponent was around 2.0, 
Therefore, the comparison of the rate constant Z could be made 
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only for the computed, values for the agitator speeds of 
300 RPM and 420 RPM, Pigixre 6»44 gives the plot of E versus 
E foi- these two agitator speeds. If the average value of 
the exponent is taken to be 1.5, 

K a 

Prom Pigure 6.44 it can be seen that the agreement 
between the computed values of K and the model is not bad 
considering the fact that there was appreciable variation in 
the value of the exponent. 

Pigure 6,45 gives the plot of computed rate constant K 
versus the interfacial tension u for region I at the agitator 
speed of 3 OO PPM. Comparison could be made only for water- 
ed^ +i-octane (50-50 por cent) and water-ihiisole + 001^(80-20 
per cent) systems for reasons stated earlier. As the value of 
the exponent n in both the cases were around 1.0, the 
functional dependence of K on a is given by, 

K a . 

The computed values compare very well with the model 

as can be seen from Pigure 6.45. Since the value of the 

exponent n for region II was around 1.8 in case of water- 

Anisole + 001^ (80-20 per cent) and water-chlorobenzene systems, 

the functional dependence of E on cr according to the model is 
—5 24 

K a a . The values of the computed rate constants K are 

more or less same in both the cases as can be seen from Table 6,10. 
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Hence, the agreement between the model and the computed 
values of K is very poor. This may be attributed to the 
peculiar behaviour of water-chlorobenzene system for which 
the region corresponding to the exponent value 1.0 could not 
be identified at small stirring times, 

6.5 EQUILIBRIUM DROP VOLUME DISTRIBUTIONS 

Equilibrium drop volxime distributions were measured for 
water-CCl^+i -octane (50-50 per cent) system at the agitator 
speeds of 420 and 300 RPM with 3 inch diameter paddle impeller 
and at the agitator speed of 480 RPM with 3' inch diameter as 
well as 2 inch diameter paddle impeller^. Measurements of 
equilibrium drop volume distributions were also made at the 
agitator speed of 300 RPM for water-Anisole + CCl^ (80-20 per 
cent) system and water-chlorobenzene system. Figures 6.46 
throu^ 6.51 give the plot of equilibrium drop volume density 
f^(d/3!^2) versus non-dimensionalized drop size d/3!^2* 

6,5.1 Prediction of the Equilibrium Drop Volume Distribution 


Solution of the population balance equation using the 
asymptotic expansion of the proposed transitional breakage 
probability under the assumption of constant coalescence 
frequency and binary -uniform breakage yielded the equilibrium 


drop volume density, 

f ( ? ) 



(3.2.19) 


whero 5 = , a = 25\ciB-±) '’ 

the density exhibits maximum. 


is the diameter for which 
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Fig. 6-46- Plot of equilibrium drop volume density fv(d/d32) 
vs. ( d/ doo ) . 
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Fig. 6-47 - Plot of equilibrium drop volume density fv(d/d32) 
VS', (d / 032 ). 
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Fig. 6-48 - Plot of equilibrium drop volume density tv (d /d32 ) 
vs. (d/dq2 ). 
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Fig. 6-50 - Plot of equilibrium drop volume density fv(d/d32) 
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Fig 6 51 - Plot of equilibrium drop volume density fv ( d / d 32 ) 
vs. ( d / d':^? ) . 
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Quantitative comparison of the derived distribution 
with the experimental data is possible only if the value of 
a is known, a is defined as 


a = (;r-) 


-5/3 


( 3 . 2 . 12 ) 


From equation (3.1.10), we get 
d 


(_S) 

'^1 ^ 


^4 


0,6 


( 3 . 1 . 10 ) 


From the correlations of Chen and Middleman [30] for 

the saute'r mean drop diameter batch agitated dispersion, 

d 

(~p) = 0.053 (6.5.1) 


i.e. 


32 

= Or- a constant. 

s 5 

It was observed by Chen and Middleman [30] that the 

d 


eqTTilibrium density exhibits a maximum at =1,06, 

^32 


1 . e. 

Therefore 


d = 1.06 d 


32 

( '.| ) = 1.06 Cp 


( 6 , . 5 . 2 ) 

(6.5.3) 


s 


Mow, the transitional breakage probability r(v) predicted 
by the model is given by 

r(v) = X erfo [{2(2^/5-1) 6^^^ -2 

^-5^/2 j 1/3] (3.6) 

Introducing the factor 0,2 in the argument of the 
complimentary error function. 
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r(v) = xerfc [0.2(2(2^^^-l) 2_)l/2 

P 

(f)l/9 V-5A8/2 5 1/3] 

r(d) = X erfc [0.2(2(2^^^-l) 

P 

d-5/6/2 e V3] 

The maximum stable drop diameter d^ can. now be defined as^ 


0.2 V’a 


2 i 17^5 

where a = 2(2^/^-l) 

r 2/5 

i,e. d-5/3 _ ^„2i.d: " ^ 


(6.5.4) 


(0.2)^(2^/^-l)6 ^ 


(6.5.5) 


Since the empirical correlation of Chen and Middleman [30] 

* 

would be used for the evaluation of ( ^ )ythe maximum stable 

s 

drop diameter d should be calculated for the vessel geometry 
and impeller used by Chen and Middleman [30]. 

Power input for a baffled stirred vessel with turbine 
impeller [15,16] is given by 

= 6.3 N^l^. 

Pg 

e = 1.27 for T=H and ^ = 0.54, the 

Jt'^h ^ 

dimensional ratios being the same as those used by Chen and 
Middleman[ 30] , Using this relation for e , we get,. 


(-f) = 0.0253 


(6.5.6) 
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Ulierefore, = 2.095 
i • e, a = 0.264 


or 



2.22 


As B = 12.25, <7^ = 23 (aCT ) "" 0*05351 (6.5.7) 

and a = 0.2313 (6,5.8) 

Therefore the equilibriiam drop volume density in terms of 


(d/a’^ 2 ) given by 

f (_!.) = ^ n 

‘^32 d* V’2Tf0.2313 


exp [~ { 


1 

1.06 


- 1 


22 


0.10702 


1 eKpf-tnoef -">^3 

f2% (0.2313) (1.06) 32 

0.10702 

(6.5.9) 

6.5.2 Comparison of the Drop Volume Density Predicted by the 
Einetic and Statistinal Thermodynamic Models : 

The equilibrium drop volume density predicted by the 
kinetic model as given by equation (6.5.9) is plotted in 
Figures 6.46 through 6.54. The equilibrium drop volume density 
predicted by the statistical thermodynamic model is also 
shown in Figures 6,46 throu^ 6,54 for comparison. It can 
be seen that the equilibrium drop volume density predicted 
by these models agree very well with each other. It is very 
interesting "feo no te that ihe equilibrium drop volume density 
arrived at through two entirely different approaches predict 
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more or less the same density. This could be viewed as 
vindication of both the models, 

^•5.3 Comparison of the Predicted Equilibrium Drop Volume 

Density with the Ezperimental Data ; 

Pigures 6.46 throu^ 6.54 compare the esperimental 
data with the predicted eq_uilibrium drop volume density. 

The equilibrium drop volume density predicted by the coales- 
cence - redispersion model [47] is also shown in Pigures 6.46 
through 6,54 for comparison. Prom Pigures 6.46 throu^ 6,51 
i-t can be seen that the experimental drop volume density 
agrees fairly well with gaussian approximation. However, 
there is a shift of the experimental data towards small non- 
dimensionalized droplet sizes d/3^2 compared to all the three 
models though the agreement with coalescence-redispersion 
model seems to be better. In Pigures 6,46 through 6,51 the 
experimental drop volume density exhibits a maximum at around 
d/d^2=0«9 in all the cases. This is in disagreement with 
the experimental observation of Chen and Middleman[ 30] who 
observed a maximum at d/3^^2 = 1.06. The maximum value of 
the density in Pigures 6.46 throu^ 6.51 is around 1,8 in all 
the cases. Figure 6,52 compares the three models with the 
experimental data of Chen and Middleman [30]. The kinetic 
and statistical thermodynamic models agree well with the 
experimental data, Pigrtre 6,53 compares the three models with 
the experimental data of Brown and Pit [36], The agreement of 
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the kinetic and statistical therraodynauiic models with the 
experimental data is good only for small sizes i.e. d/J^ 2 < 0.8. 
This is because of a lot of scatter for large drop sizes i.e. 
for d/3!j2 >0.8. The scatter can be attributed to the 
experimental errors inherent in the experimental techniq.ue 
viz, photomicrography. Photomicrography directly gives the 
number density w^ith respect to the drop diameter* Conversion 
of niimber density to volume density magnifies errors in the 
drop size measurement especially for large drop sizes as the 
conversion involves the evaluation of the third moment of the 
number density. In fact^ only Brown and Pit [ ] have 

observed a bimodel equilibrium volume density for agitated 
lean liquid-liquid dispersions. This could also be due to 
the experimental errors, Figures 6.54 compares the throe 
models with the e2;perimental data of Shinnar [53]. The agree- 
ment between the kinetic and statistical thermodynamic models 
and the experimental data seems to be good except for the 
agitator speed of 627 RPM, From PigTures 6,52 through 6,54 it can 
be observed that the coalescence-redispersion model does not 
agree with the experimental data of Chen and Middleman [30], 
Brown and Pit [36] and Shinnar [33]. Hence, it is found that 
there is a discrepancy between the present experimental 
observations of equilibrium drop volume density and the 
experimental data of Chen and Middleman [30], Brown and Pit [36] 
and Shinnar [ 53 ]. 
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The validity of this correlation for the present experimental 
data is verified. Figures 6.55 and 6,56 show plots of 
versus the agitator speed N and "versus the interfacial 
tension a. Dotted straight lines of slopes -1,2 and 0.6 are 
shown in Figures 6,55 and 6,56 respectively for comparison of 
the correlation with the present experimental data. The 
sauter mean drop diameter obtained from the present 

•flf —0 

experimental data exhibits the functional dependence d ^2 ^ 

0 6 4 

and d ^2 a cr * as can be seen from Figures 6.55 and 6,56. 
Though the values of 3^2 obtained from the present work would 
have been higher due to experimental errors, the functional 
dependence of ^ o' agreed fairly well with the 

correlation of Chen and Middleman [30], This is probably 
the reason for the consistency of the experimental equilibriim 
drop volume density f^(d/3’^2) terms of non-dimen sionalized 
drop size d/d^ 2 » 
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tor speed (RPS) 


versus agitator speed. 




CHAPTER 7 


CONCLUSIONS AND RECOMMENDATIONS 

An attempt has heen made to study drop breakage in 
agitated lean liquid -liquid dispersions. Drop breakage in 
agitated lean liquid -liquid dispersions is random and hence 
is represented by related probabilities, namely, the transi- 
tional breakage probability and the daughter droplet distri- 
bution. 

A model for transitional breakage probability of 
droplets in agitated lean liquid -liquid dispersion was 
proposed based on the mechanism of breakage of droplets due 
to their oscillations resulting from relative velocity 
fluctuations. The oscillating droplet was modelled as an 
^mdamped simple harmonic oscillator. The restoring action 
of the droplet due to the interfacial tension was represented 
by the restoring action of the spring of the simple harmonic 
oscillator, Eragmentation of the droplet was assumed to 
occur if the turbulent motion provides the minimum increase 
in the interfacial energy required for breakup. The velocity 
and pressure fluctuations occurring on the surface of a droplet 
were thought of as due to turbulent eddies * striking' the 
surface of the droplet. It was postulated that the number 
of eddies arriving on the surface of the droplet is a poisson 
process. The probability density of relative velocity between 
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two points in an agitated vessel was assumed to "be gaussian 

in tlie derivation of the transitional breakage probability, 

1 universal transitional breakage probability in terms of 

non-dimensionalized drop diameter was derived for all 

dispersed phases whose density and viscosity are almost the 

same as that of the continuous phase. The maximum stable 

drop diameter derived from the model showed a dependence 
—0 6 

of . It was shown that the proposed transitional 

breakage probability could be approximated by a series of 
'power law' models Kv^ with different exponents for different 
drop volume regions and that the value of the exponent 
increased with decrease in droplet volume. A functional 
relation for the rate constant K in the power law approxi- 
mation Kv^ of the transitional breakage probability was 
derived in terms of the parameters and physical properties 
of the system, 1 universal non-dimensional drop volume 
distribution for agitated lean liquid-liquid dispersions was 
derived by analytical solution of the population balance 
equation simplified by order of magnitude estimates under the 
assumptions of constant coalescence frequency and uniform 
binary breakage. Interestingly enou^, this analytical 
solution yielded a gaussian non-dimensionalized equilibrium 
drop volume density as suggested by Chen and Middleman [50]. 

The prediction of equilibrium drop size distribution 
by solution of population balance equation necessitates the 



knowledge of breakage and coalescence. Even for lean 
dispersions, where the interaction between droplets is 
negligible, it is necessary to take coalescence into accomit 
which makes the population balance equation non-linear. Hence 
it becomes necessary to make simplifying assumptions for 
the population balance equation to be amenable to solution. 
Prediction of equilibrium drop size distribution from a. 
statistical thermodsmamic viewpoint overcomes this drawback 
as it considers only the total energy of the system i.e, 
the kinetic, potential and iiiterfacial energies of droplets. 

A statistical thermodynamic model for equilibrium drop size 
distribution was proposed. As the interaction between droplets 
is negligible for a lean liquid -liquid dispersion, the system 
was assumed to consist of independent droplets, A canonical 
ensemble of droplets was considered in the derivation of 
equilibrium drop size distribution, Ihe energy of a droplet 
was split into”, (i) Kinetic energy of pulsating convective 
motion, (ii) Kinetic and potential energies of oscillation 
and (iii) Potential energy of a droplet at its spherical 
equilibrium shape. The energy of oscillation was calculated 
by viewing the oscillation of a droplet as forced oscillation 
of an incompressible ^here due to random periodic disturbing 
pressure acting at two diametrically opposite points on the 
surface of the sphere. Potential and kinetic energies of 
oscillation were represented in terms of a deformation 
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coordinate r] and its conjugate momentum p^. Potential 
energy of a droplet at its spherical equilihrium position 
was taken as the work done in transferring the dispersed 
phase reversibly from a bath of dispersed phase liquid at 
the same temperature, pressure and chemical potential as 
that of the continuous phase into the continuous phase in the 
form of a spherical droplet. A partition function for a 
droplet was derived. A droplet was found to have 6 degrees 
of freedom. Equilibrium drop size distribution was predicted 
from the partition function of a droplet, A universal 
non-dimen si onalized drop volume density was derived in terms 
of dimensionless drop diameter d/3^2» 

Drop volume distributions in a batch agitated lean 
liquid -liquid dispersion eventually approach an asymptotic 
distribution in the similarity variable (1+Kt)v°' if (i) the 
transitional breakage probability r(v) follows a 'power law' 
Kv^ and (ii) the breakage is 'similar', i.e. the daughter 
droplet size distribution depends only on the ratio of the 
daughter and parent droplet volumes [4,53]. The ezistence 
of SLich a similarity transformation was made use of to 
verify the model for transitional breakage probability and 
also to obtain the exponent n, the rate constant K and the 
daughter droplet distrib^ition from the transients of drop 
volume distribution in a batch agitated lean liquid -liquid 
dispersion. 
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Drop voliime distribution measurements were made in 
a jacketted glass vessel agitated by a 6 blade paddle 
impeller. Measurements were made at the dispersed phase 
volume fraction of 0.1 per cent in order to minimise drop 
coalescence. Micro encapsulation technique was employed to 
freeze the system at the time of measurement. This technique 
encapsulates the droplets with a thin nylon film formed as 
a result of an interfacial polycondensation reaction. The 
volume distribution of droplets was measured by means of 
electronic size analysis (Coulter Counter). Special probe 
assemblies consisting of platinum electrodes with aperture 
of 2 mm and 0.8 mm diameter were used to monitor the droplets. 
The pulse signals from the probe were processed and fed to 
a Multichannel analyser to give number count of droplets for 
different droplet volumes. This yielded the volume distri- 
bution of droplets. 

Transients of drop volume distribution in agitated 
lean liquid -liquid dispersions were measured for water-COl^ + 
i-Octane (50-50 per cent) system at three different agita,tor 
speeds of 480, 420 and 300 RPM with 5 inch diameter paddle 
impeller and at the agitator speed of 480 RPM with 2 inch 
diameter paddle impeller. Transients of drop volume distri- 
bution were also measured at the agitator speed of 500 RPM 
with 3 inch diameter paddle impeller for water-lnisole + CCl^ 
(80-20 per cent) and water-chlorobenzene systems. Experimental 



results showed that the transitional bi'eakage probability 
could be approximated by a 'power law' of varying exponents. 
Three distinct regions of exponents 1.0, 1.5 and 2.0 could 
be identified. The value of the exponents obtained from 
the experimental data agreed fairly well with the values 
of 1,08, 1,42 and 2,12 obtained from the approximation of 
the predicted transitional breakage probability by a series 
of 'power law' models of varying exponents. The intermediate 
value of around 1.8 obtained for the exponent for water- 
chlorobenzene and water-ini sole + ^iCl^ (80-20 per cent) 
systems in region II is believed to be due to overlapping 
of two regions. 

As the transitional breakage probability r(v) could 
be approximated by a series of 'power law' models with 
different values of exponent for different droplet volume 
range and for each of these regions Ev^*t is the similarity 
transformation which would collapse all the drop volume 
distributions in that region into a single curve, one would 
expect that a transformation r(v)‘t would collapse all the 
drop volume distributions at all times into a sin^e curve 
irrespective of their region if the model for transitional 
breakage probability is indeed valid. However, the 
computations of r(v) as predicted by the model yielded 
very large values of maximum stable drop volume v^. There- 
fore a factor of 0,2 was introduced in the arguement of the 
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comp lament ary error function in order to match the variation 
of the proposed transitional breakage probability with 
droplet volume in the range of drop volumes encountered 
with experimentally observed variation. It is believed 
that the introduction of this factor is necessary to account 
for the important inhomogeneities in the flow field of an 
agitated vessel. The transf orriiation 

erfoi (2(21/5-1) 1/5 6^/5 2 (f)l/9 v-5/lhoel/5] -t 

seems to collapse all the drop volume distributions more or 
less into a single curve thus vindicating the model for the 
transitional breakage probability. 

The drop volume distributions in a region were trans- 
formed using a modified similarity transformation 

_ V 

V and n being' an average drop volume and the exponent of 
’power law’ approximation of transitional breakage probability 
for that region respectively. The similarity transformation 
collapses all the drop volume densities in a region 

V 

into a single curve thus vindicating the assumption of 'similar 
breakage' made in the derivation of the similarity transformation. 
Hence the moments of the function Kg(x) were computed from 
the moments of the modified drop volume density f . 

V 

The function Kg(x) was then estimated from its moments through 
orthonormal polynomial expansion with the weighting function 
(3x -2x'^), The rate c onstant K of 'power law' approximation 



205 


Kv^ of the transitional breakage probability was then 

estimated from Eg(x). The estimation of the daughter droplet 

size distribution g(x) was found to be only q.ualitative which 

is believed to be due to the inclusion of insufficient number 

of terms in the orthonormal polynomial expansion. Inclusion 

of more terms in the polynomial expansion was desisted due 

to the inaccuracies involved in the evaluation of higher 

moments of modified drop volume density f { }, The 

V 

estimated daughter droplet distribution g(x) indicated that 
(i) the breakage mechanism seems to be the same for all 
systems and experimental conditions (ii) the normalised 
daughter droplet density is more or less same for all the 
cases and (iii) the daughter droplet density is broad. The 
functional dependence of the estimated rate constant K with 
respect to the agitator speed and interfacial tension agrees 
well with the model. 

Equilibrium drop volume distributions were measurod 
for water-CCl^+i -Octane (50-50 per cent) system at the 
agitator speeds of 420 and 300 RPM with 3 inch diameter 
paddle impeller and at the agitator speed of 480 EP*M with 
3 inch as well as 2 inch diameter paddle impellers. Measure- 
ments of equilibri-um drop volume distributions were also 
made at the agitator speed of 300 RPM with 3 inch diameter 
paddle impeller for water-Anisole + 001^ (80-20 per cent) 
and water— chlorobenzene systems. The measured non-dimensionalizcd 
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equilibriTjun drop volume density f-y-(d/3!^2) compared with 

the density predicted by (i) Kinetic node! (ii) statistical 
thermodynamic model and (iii) Coalescence-redispersion 
model. The non-dinensionalized equilibrium drop volume 
density pi*edicted by the kinetic and statistical 

thermodynamic models agree very well with each other ^ It 
is very interesting to note that the equilibrium drop volu'ie 
distribution arrived at through two entirely different 
approaches predict more or less the saae density. This 
could be viewed as vindication of both the models. The 
exp ex'iment al non-dim ensionalized drop volume density 
agreed fairly well with the gaussian approximation. However, 
there is a shift in the experimental drop volume density 
towards smaller dimensionless drop diameter (d/ 3 ^^ 2 ) when 
compared with the kinetic and statistical thermodynamic 
models and the agreement with the coalescence-redispersion 
model seems to be better. Comparison of .the non— dimen si onaliz 
equilibrium drop volume density predicted by these models 
with the experimental data of Chen and Middleman [30], Brown 
and Pit [56] and Shinnar [33] showed that the experimental 
data agreed fairly well with the kinetic and statistical 
thermodynamic models and that the coalescence-redispersion 
model did not agree with the experimental data. Hence it 
was found that there is a discrepancy between the present 
experimental data of equilibrium drop volimie density and the 
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experiment aJ. data of Olien and Middleman [30], Brown and Pit[36] 
and Sliinnar [55], The reason for this discrepency could be 
the loss of count of vexy small voltage pulses corresponding 
to droplets of very small volumes due to poor signal to 
noise ratio for these droplets. The vario-tion of the 
measured sauter mean drop diameter d ^2 with agitai,toi speed 
and interfacial tension agreed fairly well with the correlation 
of Chen and Middleman [50]. 

The experimental verification of the model for transi- 
tional breakage probability and the inferrence of daughter 
droplet distribution attempted in the present work is limited 
to agitated lean liquid -liquid dispersions with paddle 
impeller and a few dispersed phases. It is necessary to 
conduct similar experiments for more number of systems and 
different types of impellers in order to establish the validity 
of the present model for transitional breakage probability 
more firmly. Since the hydrodynamic characteristics of the 
flow field in a continuous well-stirred vessel is more or less 
sane as that in a batch well stirred vessel, it is believed 
that the present model for the transitional breakage probability 
should also hold good for continuous systems, fiamkrishna [4] 
has shown the method of prediction of steady state drop volTmie 
distribution in a continuous agitated lean liquid-liquid 
dispersion with the knowledge of transitional breakage 
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probability, dau^ter droplet distribution and inlet drop 
size distribution. It may be necessary to make measurements 
of steady state drop volui'ie distribution in a continuous 
agitated lean liquid-liquid dispersion with, known inlet drop 
size distribution in order to verify the validity of the 
present model for transitional breakage probability in 
continuous systems. The final objective of the study of drop 
breakage in agitated lean liquid-liquid dispersions is to 
predict rate processes in liquid -liquid dispersions by taking 
the simultaneity of heat and/or mass transfer and droplet 
phenomena into account. The present model for transitional 
breakage probability, if found satisfactory for continuous 
agitated lean liquid-liquid dispersions, could be used to 
predict the mass transfer rate of a solute from the dispersed 
phase into the continuous phase [1], in a continuous agita.ted 
loan liquid -liquid dispersion. Experiments could be conducted 
to siiiula.te a stage of a continuous multi-change extractor 
and the steady state mass transfer rates could be measured and 
compared with the theoretical predictions. 
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APPBHDIZ I 

OilDPK 01-' MAGNITUDE ESTII14TES OP TliE POPIIDATIOE BALANCE 

EQUATION 

For an equilibri-um distribution narrow around the 
maximnin stable drop diameter d^, we can divide the distri- 
bution into three zones as shown in Figure I.l. 

Zone I: Zone of droplets of anall volume, ^<< 1. 

Zone III Zone, where most of the distribution is 
concentrated,? 'v> 1. 

Zone nil Zone of droplets of large volxmie?>>l. 

For Zone I the population balance equation at equilibrium in 
non-dimensional ized volimne is 

0 = 2 f y(x') H ta' (I.l) 

X o 

as i) droplets do not break in Zone I, ii) production of 
droplets in Zone I due to coalescence of smaller drops is 
negligible. Differentiating (I.l) w.r.t. x 



For xsl, y(^) = 0 so that ^ = 0 and we write 

F(x) = nv'*x 

where n = Number of droplets in Zone I per unit interval of 
non-dimensionalized drop volume x. 

f(C) = 3nv 


Thus 


( 1 . 2 ) 
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i'roia (I.l) we have 

0 = 2 J y( C’ ) d?' 

so that dividing by dispersed phase fraction V 


CO f ( £ 

J y( C’) 

? 


,6 


d ? 


li ) T 


3 C- 


J 

o 


.3 


d C’ (1.3) 


?el 

Since f^( is concentrated around C=l 5 taking the order 
of magnitude of terms in ( 1 . 3 ) 


y( ^ ) V ^ nv* where 5 s I or II 
or 

^ £i-L xi. = ixLLl^ 

nv nv* x' 

Since nv*x' is the volnme of drops in unit non— dimensionalized 
vol-ume interval aroand x' = we may write from ( 1 . 4 ) 


or equivalently [j, » y(C )• 

Since y(? ) = A exp [- B it can easily be seen 

that for ^ s II 

Y> ( C ) 'V' Y( 5 ) '^ 20 y(C ) (1.5) 

Fow ( 1 . 6 ) 

and 


”f ( C >) 

S-^L^ as - =f -1 (1.7) 
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For Zone II, the popiilation halance equation is given by 


(3.2.7), which on division by yields 
0 


t( K )f ( 1 , , 

= — ^ „ 1 Y'(c ) f^cn 


I 


V 


II 


y(C ) 

T 

III 


+ 2 4 d 5 ' - ( 5 ) + 


IV 


V -"n 

V 


3_ILk_ nV^ 


V 

s. 

VI 

¥e now consider the coalescence integral for ? s II 


( 1 . 8 ) 


r-^ 




c5 


The contribution to this integral can come from any one of 
the following cases; 

(i) ^'s I, ^ £ IIj (ii) £ II, ? £ I; (iii) V s I 

If we denote the corresponding values of I by ^(ii)’' 

and then from (1.2) we have 


d r ^^n 
d Z 


[ 


£ (c ) 


-] 


0 


which implies that = ^(^±±±) = 0 so that I = I/',-^. 


(i) 


or 


I ^ ^ ] d r, 5’ e I, 5 e II 


, C el 
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Considering the order of magnitude of different terms in 
(1.8) it can be shown that 


I 'V 

IV'v 


iLlI. 

V » 

V 


II 

15-11 


V " 


Y 'v 


III % 

ZiLlJsSjI 5 
ay^ 15-11 


riri 

vj-v ^ Y ,. ( .. 5 . . ., I.<-5 ,.!.£ 

A 


3 

where 5=1, 5' «1, a << 1, and/^«^'''^. from these estimates 
it is clear that the coalescence integral IV is of a smaller 
order of magnitude than the other terms in (1,8). 
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APPEE'DIX II 

qalibplApiolt of the PHOBES 

Ten saiaioles of narrow distribution of spherical glass 
particles in the size range of 50 to 800 microns were used 
for the calibration of the pi'obes. Average particle volume 
of a sample was measured through microscopic counting. A 
slide of the sam;^le was prepai’ed and the mmiber counts of 
particles of different sizes were measured by viewing the 
sample through microscope. The average particle volume of 
the sample was calculated from the measured number density. 

It was necessary to calibrate the Multichannel analyser 
prior to the calibration of the probes in order to get the 
relationship between the channel number and the pulse 
amplitude. Voltage pulses of variable amplitude were fed 
from ORTEC precision pulse generator to the MCA. The pulses 
had a very small rise tiime and hence were acceptable to the 

i 

MCA. Since the voltage pulses from the precision pulse 
generator were of the same amplitude, the pulses were logged 
in one particular channel corresponding to the pulse amplitude, 

f 

Figure II. 1 shows the plot of channel number versus pulse 
amplitude. The relationship is linear as can be seen from 
Figure II, 1. By least square fit the equation of the straight 
line is 


y = 18.82 x - 0.7775 


(II. 1) 
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■where j = channel number 

and X = pulse amplitude in 'volts 

In equation (II. l) the intercept-0.7775 refers to the channel 
n'uxiiber corresponding to a pulse of zero volts. 

A sample of narro'w distribution ox spherical glass 
particles 'was suspended in a conducting medium of haCl 
solution, 'The pi'‘obe to be calibrated "was immersed in the 
suspension. The probe was filled ■with the conducting medium 
by tJithdrawing llaCl solxition into the probe through the 
application of vacuum. Current and gain settings 'vrere 
adjusted. The electrophoresis po'wer supply put on t o 
apply 60 volts D.C. bet'ween the electrodes, ho'w, a sample 
of the suspension "was withdra'vm through the aperture by slo'wly 
opening the stopcock connected to the vacuum p'ump. The 
amplitude of the amplified voltage pulses were observed in 
the oscilloscope. If the pifLses ■were not in 'bhe proper 
range, curx'ent and gain settings were readjusted to get 
voltage pulses in the proper range. The output from the 
P'ulse shaper was logged in the HCA and the spectrum of the 
sample was taken. The average channel number of the specti-um 
was calc'ulated from the printout of the MCA. Since the 
sample was a narrow distribution of glass particles, the 
average channel n-umber should correspond to the average 
particle volume of the sample. Since the pulses from the 
probe were amplified prior to the analysis by the MCA, 
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comparison of the channel members for samples of different 
average particle volumes could be done only if the amplifi- 
cation factor is maintained constant. Hence the average 
channel nrjnber of the spectrum was converted to a base 
value corresponding to a current of 1 mA and the amplifier 
gain of 100 . The same procedure was repeated for different 
samples at the same HaCl concentration to get a set of 
readings at a particular resistivity of the medium. In all 
the readings, the D.C. voltage across the electrodes at the 
reference current of 1 mA was measured. Since the resistivity 
of the medium was maintained constant for a sot of readings, 
D.C. voltage also was consta.nt. In effect, the D.C. voltage 
across the electrodes at the reference current of 1 mA 
corresponds to the resistivity of the medium. In all the 
observations, it was ensured that the average particle size 
of the sample did not exceed 40 per cent of the aperture 
diameter of the probe in order to avoid deviations from 
linear volumetric response. Since -0,7775 was the intercept 
in the calibration curve for the MCA, the set of average 
channel number and the corresponding average particle volume 
was fitted to a straight line 

y = mx - 0.7775 

where y = Average channel number 
and X = Average particle volume 
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Prom eq.uatiorL (5.2,1), ¥6 see that the amplitude of 
tlie pulse from the probe is proportional to the resistiTity 
of the medinm if the resistivity of the particle is very 
much higher than that of the medium. Therefore the same 
set of observations uere made at different haCl concentrations 
to obtain the effect of the resistivity of the medium. Since 
the I).C. voltage across the electrodes corresponds to the 
resistivit 3 '' of the medium, one ¥ould expect the slope m to 
be a linear function of the D.C. voltage, figures II. 2 and 
II. 3 give the plots of average ch8.nnel number versus particle 
volume for probes I and II for different D.C. voltages. The 
insets in figures 11,2 and 11,3 give the plot of the slope 
versus the D.C. voltage. 

The relationship between the channel number and particle 
volLimo obtained through least square fit for both the probes 
are given by. 

Probe I of aperture diameter 2 mm; 
y + 0.7775 = (0.02345 Y + 1.5041) x (II. 2) 

Probe II of aperture diameter 0.8 mm; 
y + 0.7775 = (0.16925 Y + 15.8227) x (II. 3) 

where y = Channel number 

X = Particle volume (CC x 10^) 

Y = D.C. voltage (volts) 



D.C, voltage - 16-58V 
D.C. voltage - 21-46V 



Particle volume , cc x 10° 




Fig. II. 3 -Calibration curve for probe II. 
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Since tliis relationship is obtained for the reference 
current of 1 mA, amplifier gain of 100 and unit attenuation, 
the relationship for other current, gain and attenuator 
settings could be obtained by multiplying this relationship 
by the ratio of amplification factor and the reference 
amplification factor. 
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APPENDIX III 

ElICAPSULATIOiT PECHITIQUE 

Phe encapsulation technique makes use of a polyconden- 
sa/cion reaction between a diacid chloride and a diamine. P, 
Morgan and S.L. Kwolek [64] showed that when the solution 
of a diacid chloride is brought in contact with an aqueous 
diamine solution, high polymer foi’ms in a fraction of a 
second near the interface. The reaction rates were found to 
be very high and the main reaction was found to be consi- 
derably faster than any important side reaction at the poly- 
iaerisation site. 

The chemistry of interfacial polycondensation reaction 
was found to be simple [64]. Any diacid chloride and a 
diamine are believed to react by Sjy2( nucleophilic) mechanism 
to form a protonated amide from which a proton is rapidly 
eliminated in the presence of more base. The proton acceptor 
is presumably an amine group on a diamine molecule or the 
end of an oligomer chain. Water also could act as a proton 
carrier. 

n M k 

Cl - C - (^2^3 - c'- 01 + lHGR2h ^2 “ 

0 OH 

[Cl C (CH2)g C - 1 (CH2)5 1^2 ] 

! ! 

I I 

Cl H 
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0 OH 

Cl C (CH2)g C - IT (GH2)g ^2 + Diamine salt. 

In the presence of a base (ITaOH) , the reaction occurs as 

follows ; 

0 0 

01 - C ~ (CH2)g - C - Cl + - (CH2)g - 1^2 — — 

Sebacoyl chloride Hesamethylene diamine 
^ H HO 0 

H[4 t - (CH2)g - H - C - (002)3 - 0 01 + n HaCl. 

ITylon 6--10, 

It wa.s found that the thickness of the polymer varies 
with the solvent, concentrations and reaction times. In my 
case, for very short reaction times, the thicloaess varies 
from 0.1 micron to a few microns. 

Madden and McCoy [54] used sebacoyl chloride and hexa- 
methylene diamine for encapsulation of CCl^ droplets in water, 
fhey observed that the hydrolysis of sebanoyl chloride does 
take place to some extent influencing the drop size distri- 
bution, They also observed that hydrolysis results in the 
shift of the drop size distribution towards larger size. They 
suggested that diacid chlorides of higher molecular weight 
should be used to minimise the bad effects of hydrolysis. 
Mylnek and Resnick [39 ] used terephthalic acid chloride and 
piperazine for encapsulating drops of 001^, i-octane mixture 
in water. In the present work, terephthalic acid chloride 
and piperazine were used for encapsulation. In the presence 
of a base (ITaOH), the reaction occxirs as follows. 
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Terephthalic acid chloi'-ide 


H 

I 


li 


1 


ITaOH 


pip erazine 


H 



C- 


h 


Cl + ii NaCl 
n 


The concentration of terephthalic acid chloride in 
OCl^ was 1 per cent by weight, A mixture of equal volumes of 
aqueous solutions of piperazine (2 gms/litre) and haOH 
solution (1,3^) gave diamine solution of required concen- 
tration. At the time of encapsulation, 50 cc of this 
diamine solution was added to the continuous phase. 



